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Attenuation  of  shock  waves  was  studied  in  annealed  1060  aluminum. 
2024-T351  aluminum,  and  Teflon  by  impacting  samples  with  explosively 
driven  aluminum  plates.  Free -surface  \re!ocitie«  were  measured  a*  a 
function  of  target  thickness  by  recording  the  time  of  flight  across  a  known 
distance  of  a  thin  shim  which  was  originally  in  intimate  contact  with  the 
surface  of  the  sample.  A  streak  camera  was  used  as  the  recording 
instrument.  Experimental  results  are  believed  to  be  more  accurate  than 
any  obtained  previously.  Samples  of  2024-T351  aluminum  were  shocked 
to  approximately  1 10  kbars  and  340  kbars  by  fly.  r  plates  having  velocities 
of  about  0. 125  cm/fisec  and  0.  33  cm/psec,  respectively.  1060  aluminum 
was  shocked  to  110  kbars  and  significant  differences  in  its  behavior  were 
observed  in  comparison  to  the  2024-T351  aluminum. 

Two  models  are  discussed  for  representing  elastoplastic  stress- 
strain  relations.  One  of  these  models  permits  the  yield  stress  to  be  an 
arbitrary  function  of  the  hydrostatic  pressure.  The  other  permits  both 
the  shear  modulus  and  the  yield  stress  to  vary  arbitrarily  with  the  strain. 
Results  of  calculations  using  an  artificial  viscosity  code  are  given  for 
the  two  models.  The  experimental  data  do  not  shov  a  stepwise  decrease 
of  the  free-surface  velocity  as  predicted  by  the  simple  elastic -plastic 
models  with  a  von  Mises  or  Coulomb  yield  criterion.  This  qualitative 
difference  is  attributed  to  Bauschinger  effect. 

Results  for  Teflon  indicate  that  :he  fluid  model  may  be  satisfactory 
although  the  data  are  meager  and  contain  some  inconsistencies. 


iU 


*an<m  mb  has  blub,  rmmpo m  wn  rtomt. 


t 


CONTENTS 


Section  Page 

I  INTRODUCTION  1 

II  CALCULATION  OF  SHOCK  WAVE  ATTENUATION  5 

1.  The  Constant  Poisson's  Ratio  Model  6 

2.  Variable  Shear  Modulus  Model  13 

3.  Equation  of  State  of  Explosive  Product  Gases  15 

UI  EXPERIMENTAL  METHODS  17 

1.  High  Velocity  Flyer  Plates  17 

2.  Shim  Technique  for  Free-Surface  Velocity 

Measurements  21 

3.  Fluid  Gage  25 

IV  EXPERIMENTAL  RESULTS  35 

1.,  Aluminum  (Low  Velocity  Flyer  Plates)  35 

2.  2024-T351  Aluminum  (High  Velocity  Flyer 

Plates)  39 

3.  Pressure  Dependence  of  Shear  Modulus  and 

Yield  Stress  of  2024-T351  Aluminum  39 

4.  Flow  Calculations  for  2024-T351  Aluminum  47 

5.  Teflon  49 

V  CONCLUSIONS  55 

Appendices 

I  IMMERSED  FOIL  ANALYSIS  59 

H  INDEX  OF  REFRACTION  MEASUREMENTS 

AND  DISCUSSION  OF  RELEVANT  DIELECTRIC 
THEORY  69 

REFERENCES  85 

88 


v 


DlSnilUTIOK 


ILLUSTRATIONS 


Figurc‘  g£&: 

/ 

1  Schematic  Diagram  of  Ela*toplasti( 

Stress-Strain  Relations  6 

2  Physical  Plane  for  Plate  Impact  Experiment  8 

3  Profiles  of  Stress  for  Aluminum  Projectile 

Hitting  an  Aluminum  Target  12 

4  Peak  Particle  Velocity  in  Aluminum  Target 

Hit  by  an  Aluminum  Projectile  13 

5  Cross  Section  of  Launcher  of  High  Speed 

Aluminum  Plates  and  Samples  19 

6  Trajectories  of  the  Free  Surface  of  the  High  Speed 

Flyer  Plate  and  Gas-Metal  Interfaces  20 

7  Results  of  Tests  of  Flatness  of  0.  32S  cm/psec 

Aluminum  Flying  Plates  21 

8  Plan  View  of  Attenuation  Experiment  23 

9  Typical  Record  of  Attenuation  Shot  24 

10  Water  Gage  Experimental  Arrangement  26 

11  Streak  Camera  Record  for  Shot  11, 190  28 

12  Time-Distance  Diagram  Illustrating 

the  Wave  Trajectories  in  Shot  11. 762  31 

13  Peak  Particle  Velocity  vs  Shock  Travel  (x/x  ) 

for  2024-T351  Aluminum  Impacted  at  0.  12$  cm/jjsec  35 

14  Peak  Particle  Velocity  vs  Stock  Travel  (x/x  ) 
for  2024-T351  Aluminum  and  2024  Annealed0 

Aluminum  Impacted  at  0. 125  cm/psec  37 

15  Peak  Particle  Velocity  vs  Shock  Travel  (x/x  ) 
for  1060  Annealed  Aluminum  Impacted 

at  0. 125  cm/paec  38 

16  Peak  Particle  Velocity  vs  Stock  Travel  (x/x  ) 

for  2024-T351  Aluminum  Impacted  at  0.  33  cm  'p  sec  43 

vi 


ILLUSTRATIONS  (Concluded) 


Figure 

Pag/ 

17 

Shear  Modulus,  G,  and  Yield  Stress,  (Y), 
vs  Specific  Volume,  V 

46 

18 

Comparison  of  Experimental  and  Calculated  Results, 
2024-T351  Aluminum  Impacted  at  0.  125  cm / Msec 

48 

19 

Comparison  of  Experimental  and  Calculated  Results, 
2024-T351  Aluminum  Impacted  at  0.  33  cm/psec 

48 

20 

Streak  Camera  Record,  Slot  11,864, 

Aluminum  Flying  Plate  Impacting  Teflon 

51 

21 

Calculated  and  Experimental  Results 

for  0.  31 -cm  Aluminum  Plate  Impacting  Teflon 

52 

22 

Schematic  View  of  Immersed  Foil  Gage 

59 

23 

Schematic  of  Relation  between  Shock  Front 
cud  Fluid  Cell  Components 

61 

24 

Relationship  between  Shcgck  Velocity  ^ 

Fluid  Particle  Velocity  Up.  Apparent  Velocity  q0, 

Shock  Angle  ft  and  Turning  Angle  fi 

63 

rti 


TABLES 


Table  page 

I  Value*  of  Parameter*  for  Constant  V  Stre**-Strain 

Relation*  for  Aluminum  10 

II  Equation  of  State  Parameter* 

for  Composition  B  and  Dctasheet  Explosive*  16 

ui  Plying  Plate  Attenuation  Measurements  40 

IV  Values  of  Experimental  Parameters 

for  2024-T351  Aluminum  45 

V  Molecular  Refractivitie*  of  Various  Substance* 

Calculated  from  Gaseous-  and  Liquid-State  Data  81 

VI  Results  of  Refractive  Index  Experiment*  84 

VII  Molecular  Refraction  of  Water  83 


viii 


SECTION  1 


INTRODUCTION 


An  accurate  description  of  shock  propagation  in  a  solid  requires 
knowledge  of  the  complete  constitutive  relation  among  the  stress  and 
strain  tensors  and  the  internal  energy.  Although  a  large  amount  of  data 
is  available  on  the  Hugoniot  equations  of  state  of  metals  and  other 
solids,  in  general  the  data  indicate  only  the  density  dependence  of  one 
stress  component  in  compression.  No  information  has  been  directly 
obtained  concerning  stress  wJference.  or  the  path  of  the  stress-strain 
curve  upon  release  of  pressure.  Most  available  equation  of  state  data 
are  therefore  sufficient  only  for  descriptions  a f  shock  propagation  in 
which  the  "fluid  approximation"  is  adequate. 

It  is  generally  assumed  that  at  high  pressures  material  rigidity  ia 
unimportant  so  that  a  fluid 'type  equation  of  state  suffices.  Such  a  sim¬ 
plification  seems  reasonable;  however,  little  experimental  information 
exists  to  indicate  pressure  levels  above  which  this  model  is  accurate. 

A  priori  estimates  of  &*  adequacy  of  the  fluid  mcdel  based  on  aero- 
pressure  values  of  material  strength  are  not  conclusive  because  of 
possible  increases  in  strength  with  confining  pressure. 

It  has  been  demonstrated  for  some  metals  at  lower  pressures  that 
rate -independent  elastoplastic  theory  gives  a  reasonably  accurate 
description  of  the  stress-strain  path  for  compression. 1,8  The  stress- 
strain  path  followed  on  release  of  pressure  from  s  shocked  state  is  less 
well  established  and  seems  to  be  complicated  by  the  fiauschtnger  effect.  s‘ 
Unfortuns^ly,  the  decay  of  shock  waves,  at  least  at  lower  pressures, 
may  he  sensitive  to  the  exact  chaps  of  tbs  stress -strain  curve  for  pres¬ 
sure  release,  and  direct  maths  da  fee  measuring  this  curve  are  experi¬ 
mentally  difficult. 
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The  work  reported  here  is  «?.r>  extension  of  earlier  work  in  which 
the  decay  of  shock  waves  in  several  solids  was  observed  and  compared 
with  predictions  based  on  various  cssumpi-ons  concerning  the  pressure 
release  curves.3’  6’  This  earlier  work  showed  that  rigidity  is  signifi-t 
cant  in  all  the  materials  >■  .adied  at  pressures  up  to  at  least  100  kbar. 
Agreement  between  experiment  and  theory  w  - s  obtained  in  most  cases 
by  assuming  an  elastoplastic  me  el  and  a  particular  functional  depend¬ 
ence  of  the  yield  stress  and  the  shear  modulus  on  the  pressure;  this 
representation  is  not  necessarily  unique. 

Studies  on  aluminum  during  the  reporting  period  were  aimed  toward 
obtaining  more  explicit  information  on  the  variation  of  the  shear  modu¬ 
lus  arid  the  yield  stress  at  high  pressures,  and  toward  obtaining  accurate 
knowledge  of  the  shape  of  the  pressure-release  curve  in  the  immediate 
vicinity  of  the  shocked  state.  In  this  region  the  elastic  and  plastic  re¬ 
lief  waves  are  most  distinctly  separated  and  one  can  hope  to  determine 
the  shear  modulus  and  yield  stress  independently. 

The  experiments  on  aluminum  were  conducted  on  material  in 
two  d'fferenc  li  itial  conditions.  The  major  portion  of  the  work  was  per- 
(-r  rmed  with  2024-TJ5I  in  the  as-received  condition.  Some  work  was 
also  performed  on  type  1060  to  determine  whether  the  high  pressure  be¬ 
havior  depends  significantly  on  the  initial  condition.  This  latter  mate¬ 
rial  is  aHo  of  interest  because  it  is  the  subject  of  concurrent  equation 
of  state  and  strain-rate  studies,  and  beet  use  of  Allison's  observation 
that  the  fluid  model  in  satisfactory  at  higher  pressures  in  soft 
aluminum.  8 

The  results  of  the  experiments  and  comparisons  with  the  associ¬ 
ated  theoret  c?l  calculations  on  aluminum  are  presented  in  Section  IV, 
The  model  on  which  the  calculations  are  based  is  described  in  Section  II. 
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A  few  experiments  were  conducted  on  Teflon  to  compare  its  be¬ 
havior  with  predictions  based  on  a  fluid  model.  Those  results  are 
described  in  Section  IV-5. 

Several  new  techniques  were  employed  in  these  experiments. 
Descriptions  of  the  experimental  methods  are  given  in  Section  III.  A 
method  for  propelling  reasonably  stress\free,  intact  aluminum  fiver 

v 

plates  at  a  velocity  of  0. 33  cm/fXsec  was  developed  and  used  to  produce 

\ 

controlled  shock  pulses  of  345  kbar  in  alum<num.  The  initial  free- 
surface  jump-off  velocities  were  measured  with  good  precision  by  means 
of  a  technique  that  utilizes  thin  foils  initially  in  contact  with  the  surface. 

Some  experimentation  was  performed  with  an  imrrtversed  foil  fluid  gage 

\ 

which  gave  some  additional  information  about  the  sh)s,ck  pulse  in  an 
adjacent  solid.  9  Appendix  1  shows  the  analysis  procedure  for  fluid-cell 
experiments.  In  the  course  of  analyzing  these  experiments,  some  data 
were  obtained  on  the  index  of  refraction  of  water  at  high  pressure.  These 
data  are  presented  and  discussed  with  relation  to  the  Lorenz-Lorentz 
theoretical  model  in  Appendix  II. 
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SECTION  II 

CALCULATION  OF  SHOCK  WAVE  ATTENUATION 

The  attenuation  experiments  can  be  simulated  by  the  use  of  a 
computer  code.  Two  methods  were  used  to  solve  the  flow  equations. 
One,  using  the  method  of  characteristics ,  was  restricted  to  cases  in 
which  rigidity  was  neglected.  The  equation  of  state  used  in  this  code  is 


where  P  is  the  pressure,  p  is  the  density,  and  is  the  density  at 
zero  pressure.  Values  of  the  eomtants  A  and  y  for  aluminum  are 
0.  196  Mbar  and  4.  1,  respectively,  and  the  initial  density  is  2.  785  g/cc. 
Results  of  the  characteristics  code  and  Eq.  (1)  are  labeled  "fluid"  in 
the  figures,  to  emphasize  that  rigidity  was  neglected  in  the  calculations. 

The  other  method  for  solving  the  flow  equations  made  use  of  the 
method  developed  by  von  Neumann  and  Richtmyer  which  uses  an  arti¬ 
ficial  viscosity  to  smooth  discontinuities  in  the  flow.10  This  method  is 
more  easily  applied  to  problems  in  which  the  flow  must  be  calculated 
across  interfaces  between  different  materials  than  is  the  method  of 
characteristics.  It  is  also  more  easily  applied  to  problems  involving 
more  complicated  equations  of  state,  such  as  an  elastoplastic  relation. 
Situations  where  rigidity  is  neglected  are  handled  by  the  use  of 
the  equation 

P  =  Ap  +  bp  +  Cpa  (2) 

where  p  =  plp0  -  I-  In  such  a  case,  the  Hugoniot  and  the  expansion 
adiabats  are  assumed  to  coincide  and  are  all  described  by  Eq.  (2). 
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1. 


The  Constant  Poisson's  Ratio  Model 


The  elastoplast ic  relation  is  diagrammed  in  Fig.  1.  The  hydro¬ 
static  curve  is  represented  by  Eq.  {2),  and  the  upper  and  lower  curves 
are  given  by 

ax  =  P  ±  |  Y  (3) 

where  &x  is  the  stress  in  the  direction  o£  propagation  of  the  shock  and 
Y  is  the  yield  stress  in  simple  tension,  i.e.  ,  twice  the  maximum  re¬ 
solved  shear  stress.  The  upper,  or  loading,  curve  is  made  to  coincide 
with  the  Hugoniot  curve  just  at  Eq.  (1)  was  forced  to  do.  Calculated 
results  agreed  more  closely  with  experimental  results  when  Y  was  made 
to  vary  with  the  hydrostatic  pressure  as 


FIG.  1  SCHEMATIC  DIAGRAM  OF  ELASTOPLASTIC  STRESS- 
STRAIN  RELATIONS 
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Y  =  Y  +  M(P-P, ) 
o  '  a' 


(4) 


where  Y^  is  the  initial  yield  stress,  P&  is  as  defined  in  Fig.  1,  and  M 
is  a  constant.  Values  of  the  various  parameters  for  the  elastoplastic 
equation  of  state  are  given  in  Table  I  for  aluminum. 

In  the  derivation  of  the  elastoplastic  relations,  stress  0X  is  re¬ 
lated  to  density,  p,  by 


for  an  elastic  event.  In  Eq.  (5),  K  is  the  bulk  modulus,  G  is  the 
rigidity  modulus,  p  is  the  density,  and  the  subscript  x  has  been  dropped. 
In  the  constant  Poisson's  ratio  model,  K  is  replaced  by  -  V  dP/dV,  where 
V  is  the  specific  volume,  P  is  the  hydrostatic  pressure,  and  G  is  re¬ 
placed  by 


(l-2t/)  _  3V(l-2y)  dP 

2(1+*/)  "  '  2(1+1/)  dV 


(6) 


where  V  is  Poisson's  ratio.  (This  model  was  used  in  most  of  the  cal¬ 
culations  made  earlier  in  this  and  preceding  projects.  )  Combining 
Eqs.  (5)  and  (6)  gives 


do  _  a  _  dP 

dp  _ (l+*4  dp 


where  c  is  the  sound  speed  and  dP/dp  is  the  slope  of  the  hydrostat. 
Using  Eq.  (3)  then  gives 


il  m* i 

[“  do  2 

— 1 

¥  u+v) 

[dp  *  3 

dp  J 

so  that  the  elastic  sound  speed  depends  on  both  v  and  Y.  This  is  the 
speed  of  sound  in  the  shocked  material  and  is  associated  with  the  head 
of  the  rarefaction  wave  BM  shown  in  Fig  2,  The  actual  velocity  of 
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FIG.  2  PHYSICAL  PLANE  FOR  PLATE  IMPACT  EXPERIMENT 


this  wave  is  u+  c  ,  where  u  is  the  particle  velocity  in  the  shocked 
region.  For  flyer  plate  experiments,  the  particle  velocity  is  deter¬ 
mined  either  by  measuring  the  free-surface  velocity  of  a  thin  target,  or 
by  measuring  the  flyer  plate  velocity.  The  latter  is  more  desirable  in 
principle  because  u  is  exactly  one-half  the  flyer  plate  velocity,  while 
it  is  only  approximately  half  the  free-surface  velocity. 11  Experiment¬ 
ally  it  is  easier  to  measure  the  free-surface  velocity.  It  can  be  shown 
that  the  sound  speed  is  determined  if  the  physical  coordinate  of  the 


point  M  (Fig.  2)  is  known, 


where  d  is  the  thickness  of  the  flyer  plate,  x  is  the  physical  coordinate 
of  the  point  M,  and  U  is  the  velocity  of  the  shock  front.  Because  the 
particle  velocity,  u,  is  known,  the  shock  velocity,  U,  is  determined 
from  the  Hugoniot  relations.  Equations  (8)  and  (9)  are  important  links 
between  experimental  observations  and  the  theory. 


The  quantity  d&/dp  in  Eq.  (8)  is  the  slope  of  the  Hugoniot  curve 
from  which  the  sound  speed  was  obtained  when  calculations  were  done 
with  the  method  of  characteristics.  That  is,  it  was  assumed  that  the 
expansion  adiabats  coincided  with  the  Hugoniot  curve  for  those 
calculations.  This  assumption  is  also  made  in  most  impedance  mis¬ 
match  calculations,  which  are  frequently  used  in  obtaining  data  for  the 
Hugoniot.  It  is  of  interest  to  compare  the  elastic  wave  sound  velocity. 
ce  ,  as  given  by  Eq.  (8)  and  the  sound  velocity  from  the  Hugoniot.  . 
The  ratio  is 


(10) 


where  dY /do  should  be  evaluated  at  a  point  such  as  e  on  the  upper 
curve  in  Fig.  i.  This  means  that  c  is  the  speed  of  sound  at  the  head 
of  the  elastic  wave,  say  the  line  BM  in  Fig.  2.  The  derivative  can  be 
evaluated  by  using  Eqs.  (3)  and  (4);  the  result  is  M/{  1+2  M/ J).  For 
the  values  of  M  given  in  Table  1,  the  derivative  is  small  compared  to 
1.0.  Equation  (10)  then  gives  a  ratio  of  about  1.22  when  V  -  1/3.  and 
the  ratio  is  independent  of  the  shock  strength. 

Integration  of  Eq.  (S)  with  the  assumptions  given  above  gives 
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Table  I 


VALDES  Or  RUUMXTIBS  FOR  COKSTAJH 
STRESS  -STBAIE  RHATI0K8  FOR  ALDHII 


Parameters 

Variable 
Yield  Model 

Cons tan t 
Yield  Model 

Y  (Mbars) 

0.0025 

0.0025 

M 

0.055 

0.0 

P  (g/cc) 

2. TBS 

2.7B5 

A  (Mbars) 

0.755 

0.743 

B  (liters) 

1.29 

1.74 

C  (Mbars) 

1.197 

0.329 

ae  '  °f  = 


.  .llr.Ji.l  (y  +  y  ) 
(i  -  21/)  '  e  V 


dir 


where  the  subscripts  refer  f.o  the  points  e  and  f  in  Fig.  1  and  Y  is  thi* 
yield  stress.  For  V  -  i/3,  and  if  Y  does  not  depend  strongly  on  the 
strain,  the  usual  result 


a{  -  4Y 


(1 D 


is  obtained. 


The  model  described  above  has  been  used  to  calculate  the  stresses 
induced  *nen  an  aluminum  plate  hits  an  aluminum  target.  When  the 
velocity  of  the  projectile  is  0.  125  cm/psec,  the  induced  stress  is  about 
110  kbar  and  the  calculations  compare  favorably  with  experimental 
results.  Hence  there  is  some  validity  to  the  model.  There  are,  how  • 
ever,  some  objections  to  it.  One  is  the  fact  that  the  yield,  Y,  continu¬ 
ally  increases  with  the  hydrostatic  pressure.  With  a  constant  value  of 
V,  the  elastic  wave  amplitude  then  continues  to  increase,  so  that  fluid 
behavior  is  never  approached.  Another  objection  is  the  fact  that  C.  the 
shear  modulus,  also  continues  to  increase  as  -  V  dP/dV  or  pdP/dp,  so 
that  there  is  no  approach  to  fluid  behavior. 

Some  results  typical  of  those  obtained  with  the  constant  v  model 
and  the  Q-eode  are  given  in  Fig.  3.  These  results  are  for  the  case  of 
an  aluminum  projectile  0.  125 -inch  thick  hitting  a  semi-infinite  target. 

The  pressure  versus  distance  profiles  are  given  at  intervals  of  1/2  psec 
following  projectile  impact.  The  parameters  given  in  column  1  of 
Table  1  were  used  in  the  elastoplastic  stress-strain  relations.  In 
Fig.  3,  the  elastic  relief  wave  reduces  the  amplitude  of  the  pressure 
wave  by  about  30  kbar.  Similar  profiles  of  the  particle  velocity  may 
be  obtained  from  the  calculations.  Figure  4  shows  only  the  envelope 
of  such  particle  velocity  profiles,  along  with  the  results  of  the  charac¬ 
teristic  code  used  with  Eq.  (1).  Comparison  of  the  two  sets  oi  result* 
shows  the  early  attenuation  which  results  when  the  elastoplastic  stress  - 
strain  relations  are  used.  Experimental  results  from  two  previously 

*  y 

reported  experiments  are  included  »n  the  figure. 
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FIG.  4  PEAK  PARTICLE  VELOCITY  IN  ALUMINUM 

TARGET  HIT  BY  AN  ALUMINUM  PROJECTILE 


2.  Variable  Shear  Modulus  Model 

It  is  not  necessary  to  keep  the  value  of  Poisson'b  ratio  constant. 

If  V  is  permitted  to  increase  with  stress,  G  ,  the  shear  modulus  changes 
with  stress  in  a  different  way  from  that  used  in  the  earlier  calculations 
(see  Eq.  [6]).  Better  understanding  of  the  relations  may  result  if  U  is 
eliminated  and  G  is  retained.  Hence  from  Eq.  (5)  the  elastic  sound 
speed  is 


where  F  is  called  the  longitudinal  elastic  modulus.  Experiments  with 
flyer  plates  give  values  of  both  t  and  V,  so  that  F  may  be  calculated. 
If  experimental  results  are  available  at  two  or  more  stress  levels,  the 
behavior  of  F  with  stress,  or  w.th  strain,  may  be  ascertained.  It  can 
then  be  assumed  that  K  s  - V dP/dV  as  before,  so  that 

G  -  -i(F-K»  =  -A  "«) 

The  quantity  dP/dV  is  not  known,  but  because  the  stress,  a  .  >s  knu*n 
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in  each  set  of  experiments,  it  can  oe  approximated  by  dCpj/dV  where 
Opj  is  on  the  upper,  cr  Hugoniot,  curve  of  Fig.  1.  In  this  way  the  vari¬ 
ables  G.  F,  and  K  may  be  evaluated  from  knowledge  of  experimentally 
related  values  of  p  and  c. 


Comparison  of  the  elastic  and  hydrodynamic  sound  speeds  is  done 

2 

as  before  by  replacing  the  hydrostatic  pressure  with  \  so  that 

Eq.  (13)  becomes 


2  _dY_  4  _G 

dp  ’  3  dp  3  0 


(15) 


Hence  the  ratio  of  sound  speeds  is 

Ce  .  1 1  t  (a  G  2  «\!<  (16) 

ch  r  V  P  '  3  *>)/  ** 


where  the  term  involving  G  dominates  the  term  containing  Y  for  low 
values  of  the  stress.  The  results  given  in  Section  IV  indicate  that  the 
ratio  is  about  1.22  for  a  110-kbar  shock  and  about  1.  16  for  a  34C-kbar 
shock.  This  variable  ratio  is  a  desirable  feature  of  the  new  model.  The 
yield  and  the  shear  modulus  could  be  made  to  vary  in  such  a  manner 
that  the  ratio  becomes  1.  0  at  some  stress.  For  the  present,  tempera' 
ture  is  not  included  explicitly.  Effects  of  temperature  are  implicit  in 
the  functional  relationships  of  Y  and  G  to  the  density. 

♦  , 

Equation  (11)  now  becomes 


V  *  (Ye+Yf)  ^k4g ^2G 


(17) 


or  by  the  use  of  Eq.  (13). 

Y 


■V  Y<  •  '  H1— j 


(18) 


where  K  can  be  approximated  as  explained  above.  It  was  expected  that 


14 


the  experiments  which  give  p  and  c  for  an  elastic  wave  would  also 
give,  at  least  approximately,  values  of  (<ye -CT^) >  so  that  the  value  of 
(Ye+Yf)  could  be  calculated.  Once  these  values  are  known  aa,  say. 
functions  of  the  volume,  Eq.  (3)  can  be  used  to  construct  a  tentative 
hydrostat,  and  the  process  of  calculating  G  and  ( Y e+ Y £>  can  bo 
repeated.  Results  of  such  a  process  are  given  in  Section  IV.  In  antici¬ 
pation  of  the  discussion  of  the  results  it  can  be  stated  that  the  shear 
modulus  does  not  appear  to  increase  indefinitely,  a  desirable  charac¬ 
teristic  of  the  model.  Because  the  experiments  fail  to  show  a  definite 
separation  of  the  elastic  and  plastic  relief  waves,  the  drop  in  stress 
caused  by  the  elastic  relief  wave  is  poorly  determined.  This  means 
that  the  variation  of  the  yield  stress  with  stress  or  strain  is  also  poorly 
determined. 

Results  obtained  with  the  Q-code  when  the  variable  shear  modu¬ 
lus  mode',  was  used  are  given  in  Section  IV.  Comparison  with  experi¬ 
mental  results  shows  that  this  model  is  also  incapable  of  describing 
the  release  of  stress  properly.  The  variation  of  both  the  shear  modulus 
and  the  yield  stress  is  also  discussed  further  in  Section  IV. 

3 .  Equation  of  State  of  Explosive  Product  Gases 

The  Q-code  permits  the  calculation  of  the  flow  in  a  detonating 

1  3 

explosive  by  use  of  the  burn  fraction  described  by  Wilkins.  For  the 
calculations  used  here  in  the  development  of  the  high  velocity  flyer 
plate,  a  polytropic  equation  of  state  was  used  lor  the  detonation  product 
gases: 

P  =  (y-l)E/V  (19) 

where  E  is  the  internal  energy  and  V  is  the  specific  volume.  Values 
of  the  parameters  for  both  Composition  B  and  for  Du  Pont  sheet  explos¬ 
ive  EL-506D  are  given  in  Table  II. 
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SECTION  III 


EXPERIMENTAL  METHODS 

Several  new  techniques  were  developed  and  used  in  accomplishing 
the  goals  of  the  program.  These  include  a  method  for  accelerating 
aluminum  plates  to  a  velocity  of  0.33  mm/fisec  with  an  explosive  arrange  ¬ 
ment' that  permits  the  plate  to  be  reasonably  stress-free  at  the  time  of 
impact  on  the  target.  In  addition,  a  recording  system  was  developed 
that  permits  more  accurate  measurement  of  the  initial  jump-off  velocity 
of  a  free  surface  than  could  be  achieved  previously.  Finally,  attempts 
were  made  to  adapt  a  fluid-cell  gage  developed  on  another  project9  to 
the  measurement  of  particle-velocity  profiles  in  aluminum;  this  effort 
was  not  successful.  Detailed  descriptions  of  each  of  these  methods 
are  presented  below. 

1 .  High  Velocity  Flyer  Plates 

Most  of  the  previous  Stanford  Research  Institute  work  on  attenu¬ 
ation  of  shock  waves  has  used  aluminum  plates  whose  velocities  were 
about  0.  12  to  0.  13  cm/fisec.  Impact  by  these  plates  produces  stresses 
of  about  110  kbars  in  aluminum  targets.  Some  plates  weie  also  thrown  at 
0.  19  cm/(isec;B  later  it  was  found  that  these  plates  were  spalling.6 
One  of  the  tasks  in  the  current  effort  was  to  find  a  means  of  projecting 
aluminum  plates  at  higher  velocity  in  a  reliable  manner.  Such  plates 
would  permit  experimentation  at  higher  stress  levels  so  that,  for  ex¬ 
ample,  the  behavior  of  the  elastic  moduli  with  pressure  could  be  deduced. 

The  0.  19  cm/fisec  plates  spalled  partly  because  of  the  impedance 
mismatch  between  the  aluminum  and  the  brass  plate.  This  mismatch 
had  the  useful  effect  that  the  aluminum  separated  from  the  brass  to 
produce  a  free  flying  plate.  Balchan  and  Cowan14  had  thrown  plates 
from  a  high  impedance  driver  plate  by  putting  a  thin  layer  of  a  low 
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impedance  material,  such  as  magnesium,  between  the  two  plates.  This 
method  was  tried  using  a  2-inch-thick  charge  of  Comp  B  to  drive  a 
sandwich  of  0.  5-inch-thick  brass,  0.  125-inch-thick  magnesium  (or 
Plexiglas),  and  the  0.  125-inch~thick  aluminum  plate.  The  plates  were 
tested  for  spalling  by  causing  them  to  impact  Manganin  wire  gages.  The 
records  showed  the  flat-topped  wave  profile  characteristic  of  plate 
impact.  However,  the  pulse  was  not  long  enough  in  time,  and  a  second 
pulse  closely  followed  the  first.  This  was  interpreted  to  mean  that  the 
plates  had  spalled.  When  the  magnesium  or  Plexiglas  layers  were  re¬ 
duced  to  about  1/ 16-inch  thick,  spalling  was  again  observed. 

It  was  proposed  that  the  lew  impedance  layer  between  the  driver 
and  the  flyer  plates  be  replaced  by  explosive.  This  situation  was  inves¬ 
tigated  by  the  use  of  a  computer  code  using  the  artificial  viscosity 
method.  It  was  soon  apparent  that  the  brass  plate  would  have  to  be  re¬ 
placed  with  a  softer  material.  When  aluminum  was  used  as  the  driver 
plate,  it  moved  into  the  thin  layer  of  explosive  at  a  velocity  close  to  the 
Chapman -Jouguet  (C-J)  particle  velocity.  The  Q-code  indicated  that  the 
aluminum  projectile  plate  might  spall,  apparently  due  to  the  release  of 
pressure  from  the  back  of  the  explosive  charge- -the  Taylor  wave.  The 
character  of  the  relief  from  the  back  of  the  charge  was  changed  rather 
drastically  by  spacing  the  2-inch-thick  charge  of  Comp  B  1.0  cm  from 
the  0.  5-inch-thick  aluminum  driver  plate  (see  Fig.  5).  The  code  indi¬ 
cated  that  this  arrangement  produced  a  flyer  plate  that  remained  intact, 
and  that  it  had  a  velocity  cf  0.  36  cm/|isec.  The  assembly  was  checked 
by  doing  an  experiment  using  Du  Pont  Detasheet  EL-506D,  0.  125 -inch 
thick.  The  Manganin  wire  gage  records  showed  that  the  plates  were 
not  spalling  and  that  the  velocity  of  the  plates  was  0.  32  cm/fisec. 

In  the  attenuation  studies  the  plates  were  allowed  to  move  about 
1.0  inch  before  impacting  the  targets.  This  permitted  waves  in  the 
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FIG.  5  CROSS  SECTION  OF  LAUNCHER  OF  HIGH  SPEED  ALUMINUM  PLATES 
AND  SAMPLES 

plate  to  reverberate  about  eight  times,  so  that  approximately  stress- 
free  conditions  were  reached.  The  explosive  gases  continued  to  exert 
some  pressure  on  the  back  of  the  flyer  plates.  The  acceleration  due 
to  this  residual  pressure  has  not  been  measured.  It  might  have  been 
desirable  to  let  the  plates  move  farther  so  that  waves  could  attenuate 
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farther  and  so  that  the  gas  pressure  would  decrease  farther.  However, 
the  plates  were  not  flat  and  might  have  been  changing  shape.  The 
l. 0-inch  travel  was  therefore  a  compromise. 

Figure  6  shows  some  of  the  results  of  the  calculations  for  the 
flyer  plate  takeoff  described  above.  In  the  figure,  the  positions  of  the 


0«T4*Cf  FROM  VOi  OF  COMP  S  CHMMC — cm 

FIG,  6  TRAJECTORIES  OF  THE  FREE  SURFACE  OF  THE  HIGH  SPEED  FLYER 
PLATE  AND  GAS-METAL  INTERFACES  {Results  or.  from  Q-coda  calculation. 
Main  charge  was  5-<m-»hiclc  Comp.  B.) 


surfaces  of  the  two  aluminum  plates  are  given  as  functions  of  time. 

The  flyer  plate  velocity  changes  very  little  after  the  first  microsecond. 
The  shapes  of  two  plates  are  shown  in  Fig.  7,  The  plates  were  origin¬ 
ally  6  inches  in  diameter;  their  shape  is  inferred  from  a  streak  camera 
record  of  the  impact  of  the  plates  with  flat  glass  targets.  For  one  test, 
the  center  of  the  plate  arrived  at  the  glass  witness  plate  0.04  psec 
ahead  of  that  part  of  the  plate  on  a  4 -inch  diameter.  For  another  test. 
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FIG.  7  RESULTS  OF  TESTS  OF  FLATNESS  OF  0.325  em/jusc 
ALUMINUM  FLYING  PLATES 

/ 

the  difference  in  arrival  time  was  0.06  tiaec.  These  differences  in 
time  correspond  to  a  deformation  of  about  0.005  to  0.008  inch,  respec¬ 
tively.  across  a  4-inch  plate. 

2.  Shim  Technique  for  Free-Surface  Velocity  Measurements 

The  free- surface  velocity  of  a  shock-loaded  sample  can  be  deter¬ 
mined  by  several  different  methods.  1#  One  of  the  simplest  is  to  record 
the  time  of  flight  of  the  free  surface  across  a  gap.  This  method  has  the 
disadvantage  that  it  gives  an  average  velocity  in  those  cases  in  which 
the  shock  is  not  a  uniform  shock,  i.e. ,  the  pressure  profile  is  not 
flat-topped.  What  is  wanted  from  the  measurement  is  the  velocity  of 
the  surface  at  the  instant  of  reflection  of  the  shock  wave.  The 
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measurement  can  be  made  more  accurately  if  a  thin  shim  is  held  in 
contact  with  the  specimen.  Because  the  shim  is  made  of  the  same 
material  as  the  specimen,  it  acquires  the  same  velocity  as  its  free 
surface.  If  there  is  any  attenuation  of  the  shock,  the  surface  of  the 
specimen  is  decelerated  while  the  shim  continues  at  uniform  velocity. 
Four  samples  of  different  thicknesses  were  used  in  each  shot,  as  shown 
in  the  plan  view  of  Fig.  8.  The  figure  shows  the  area  of  each  sample 
that  was  viewed  by  the  camera  through  the  slit  and  the  reflections  of 
lines  in  the  mask  that  covers  the  light  source.  A  cross  section  of  the 
arrangement  of  die  charge  and  samples  is  shown  in  Fig.  5. 

For  aluminum  and  copper,  the  reflectivity  of  the  shim  changes 
sufficiently  when  it  is  accelerated  by  the  shock  so  that  the  take-off  can 
be  observed  on  the  smear  camera  record  (see  lines  A -A  in  Fig.  9). 

The  gap  is  defined  by  a  mirror  (which  may  be  partially  coated  with  gold) 
set  at  a  known  distance  from  the  original  position  of  the  shim.  When 
the  shim  collides  with  the  mirror,  another  change  of  reflectivity  occurs, 
so  the  arrival  can  be  observed  in  die  record  (lines  B-B  in  Fig.  9).  The 
elastic  precursor  wave  in  as-received  2024-T351  aluminum  does  not 
change  the  reflectivity  of  the  shim  sufficiently  for  the  precursor  to  be 
observed.  The  gap  closure  caused  by  the  unobserved  elastic  wave 
causes  negligible  error  in  the  measurements  When  annealed  2024-T351 
aluminum  specimens  whose  thickness  was  0.  125  inch  were  hit  by  a  flyer 
plate  of  the  tame  material  (unannealed),  the  velocity  as  recorded  by 
the  shims  was  essentially  the  same  as  the  velocity  of  the  projectile 
plate.  In  this  case  the  gap  was  about  0. 19  cm,  so  the  shim  was  in 
motion  for  about  1.5  psec.  During  this  time,  relief  waves  from  the 
back  of  the  projectile  plate  overtook  the  front  surface  of  the  specimen 
and  slowed  it,  The  record  (Fig.  9)  shows  the  arrival  of  the  shims  at  the 
mirrors,  lines  B-B,  and  a  fraction  of  a  microsecond  later  the  arrival 
of  the  surface  of  the  target,  lines  C-C.  This  illustrates  the  utility  of 
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FIG.  8  PLAN  VIEW  OP  ATTENUATION  EXPERIMENT  (Four  somples  are  included  in  each  *hn?t 


the  use  of  shims  ior  the  measurement  of  the  free-surface  velocity  of 
aluminum  specimens.  The  delay  between  the  arrival  of  the  shim  and 
the  target  surface  at  the  mirror  depends  on  the  thickness  of  the  target, 
i.e. ,  on  the  amount  the  shock  wave  has  attenuated. 

In  the  pa*i,  targets  in  the  shape  of  wedges  were  used  so  that  the 
free-surface  velocity  could  be  obtained  in  one  experiment  for  a  range 
of  distances  of  shock  travel.  The  optical  lever  arm  technique  was  used 
to  record  the  motion  of  die  free  surface.  One  criticism  of  the  technique 
is  that  the  shock  is  obliquely  incident  on  the  free  surface  so  that  the  re¬ 
flection  of  theteock  results  in  shear  waves.  The  analysis  is  aiso  some¬ 
what  complicated,  requiring  the  differentiation  of  numerical  data.  If 
truly  flat  flyer  plates  could  be  produced,  differentiation  of  the  data  would 
be  greatly  simplified  and  analysis  of  the  streak  ramera  records  would 
be  straightforward.  Okie  ef  che  reasons  for  the  use  of  die  snim  technique 
in  this  program  was  to  check  results  obtained  previously  with  the  optical 
lever  arm  technique. 

3.  Fluid  Gage 

a.  Technique 

Shock  attenuation  studies  are  concerned  with  die  pressure 
or  particle  velocity  profile  of  the  shock  waye  as  it  changes  with  time  in 
various  materials.  Previously  the  profile  has  been  inferred  with  the 
use  of  both  electrical  gapes  and  measurements  of  free-surface  velocity. 
Use  of  the  gages  is  desirable  ia  general  because,  unlike  free-surface 
measurements,  they  five  information  about  the  material  while  the 
pressure  ir  the  material  it  substantially  above  atmospheric  The  de¬ 
velopment  of  the  immersed  Mylar  foil  gage  provides  in  principle  a 
means  of  measuring  the  shock  profile  in  the  liquid  of  the  gage  and  also 
to  some  extent  the  profile  in  the  driving  substance . 
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The  gage  is  constructed  by  cementing  an  aluminized  Mylar  foil 
(0.09025-  to  0.0005-inch  thick)  between  rings  of  lucite  (see  Fig.  10). 

This  cell  is  mounted  on  the  material  dirough  which  the  shock  will  prop¬ 
agate  and  is  oriented  so  that  the  reflecting  foil  is  inclined  to  the  speci¬ 
men  surface  and  hence  the  shock  front  at  an  angle  Of.  The  cell  is  covered 
aith  a  thin  glass  cover  and  the  fluid  is  added  so  that  the  foil  is  com¬ 
pletely  immersed. 
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FIG.  10  WATER  GAGE  EXPF***EHTAL  AftRANBEMENT 


The  fluid  motion  behind  the  shock  front  imparts  a  velocity  to  the 
foil  parallel  to  die  shock  front  velocity.  As  a  result,  the  foil  is  bent 
through  an  angle  fi  toward  the  shock  front.  The  relationship  cf  this 
turning  angle  tc  the  fluid  velocity  is  discussed  in  Appendix  I. 
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In  practice,  the  optic  axis  is  aligned  perpendicular  to  the  unturned 
loi'  Small  light  sources  are  provided  by  placing  a  grid  of  alternately 
opaque  and  transparent  lines  over  an  extended  light  source.  The  camera 
photographs  the  images  of  these  grid  lines  in  the  highly  reflecting  foil. 

Figure  11  is  an  example  of  the  streak  camera  record  obtained  us¬ 
ing  this  technique.  As  the  shock  moves  into  the  reflecting  foil,  the  foil 
bends  and  the  grid  sources  are  then  imaged  in  both  the  stationary  and 
the  moving  foil.  This  results  in  a  jump  of  the  trace  on  the  film.  The 
magnitude  of  the  jump  and  its  relationship  to  the  particle  velocity  in 
the  fluid  are  derived  in  Appendix  I.  The  position  of  the  displaced  trace 
is  not  only  affected  by  the  angle  /?  through  which  the  foil  has  turned  but 
also  by  the  refractive  index  of  the  shocked  fluid.  This  is  because  the 
light  rays  are  refracted  as  they  pass  twice  through  the  shock  front.  Also 
of  importance  is  the  velocity  of  the  point  of  intersection  between  the  foil 
and  the  shock  front.  Its  velocity  depends  on  the  shock  velocity  and  the 
foil  angle  and  is  indicated  by  the  slanted  line,  AB,  across  the  film 
record. 

The  immersed  foil  gage  can  be  utilized  for  several  measurements, 
among  which  are:  peak  shock  pressure  in  the  liquid;  peak  shock  pres¬ 
sure  in  the  driver;  determination  of  the  refractive  index  of  the  shocked 
liquid;  and  observation  of  overtaking  shocks  and  relief  waves.  It  may 
also  be  possible  to  interpret  the  motion  of  the  grid  lines  after  jump  in 
order  to  obtain  information  about  the  shock  profile  as  a  function  of  time. 

The  peak  shock  pressure  in  the  liquid  can  be  determined  provided 
the  refractive  index  is  known.  Since  the  foil  is  inclined,  the  pressure 
can  be  measured  as  a  function  of  depth  in  the  liquid.  The  reduction  of 
the  pertinent  data  is  described  in  Appendix  I. 
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The  peak  shock  pressure  in  the  driver  can  be  inferred  with  the 
use  of  the  ceil  if  the  wave  is  not  attenuating.  The  peak  pressure  in  the 
fluid  is  calculated  using  measurements  obtained  in  ^he  fluid,  and  the 
corresponding  peak  shock  pressure  in  the  driver  is  determined  by  im¬ 
pedance  matching.  However,  if  the  shock  wave  is  not  flat-topped,  the 
pressure  measured  at  the  foil  will  in  general  not  be  representative  of 
the  peak  pressure  in  the  driver,  since  the  peak  pressure  changes  as  the 
front  propagates  through  the  liquid. 

Analysis  of  the  inclined  foil  record  gives  the  pressure  at  the 
several  locations  of  the  jumped  grid  lines.  Since  the  foil  is  inclined 
with  respect  to  the  direction  of  shock  propagation,  the  peak  shock  pres¬ 
sure  as  a  function  of  depth  U  obtained.  With  the  use  of  large  foils  or 
double  foils  (i.  e. ,  two  foils  per  cell)  the  peak  pressure  can  be  measured 
over  considerable  depth  in  the  liquid.  This  information  is  relevant  to 
the  attenuation  characteristics  of  the  fluid. 

In  Fig.  11  it  can  be  noted  that  the  grid  traces  are  visible  after 
they  have  jumped.  Since  the  foil  moves  with  the  fluid  velocity,  the  mo¬ 
tion  of  the  displaced  grid  lines  should  indicate  the  motion  of  the  fluid 
after  the  shock  front  has  passed  over  the  foil.  In  principle,  the  gage  is 
then  capable  of  providing  information  not  only  about  the  peak  shock  pres¬ 
sure  but  also  about  a  portion  of  the  relief.  This  type  of  information-- 
l.a.  ,  pulse  shape  in  the  liquid  as  a  function  of  time  and  depth- -can  then 
serve  as  a  check  of  computer  calculations  of  the  shock  flow  in  the  fluid 
and  in  the  metal  cell  bottom. 

The  refractive  index  is  known  to  depend  on  the  density  and,  less 
importantly,  on  the  temperature  of  a  liquid.  The  gage  can  be  employed 
to  determine  the  refractive  index  of  shocked  liquids .  In  fact,  a  knowl¬ 
edge  of  the  refractive,  index  is  necessary  before  any  of  the  shock  flow 
parameters  can  be  obtained  with  this  device. 
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The  liquid  cell  can  also  be  designed  so  that  disturbances  overtak¬ 
ing  the  shock  from  the  rear  can  be  observed  with  the  inclined  foil 
arrangement.  The  propagation  of  shocks  or  relief  waves  is  indicated 
by  additional  motions  of  the  displaced  traces.  In  the  case  of  overtaking 
shocks,  the  grid  lines  jump  a  second  time  as  the  shock  passes  over  the 
foil,  and  the  motion  of  the  second  shock  front  is  indicated  by  another 
slant  line  on  the  film. 

Figure  11  is  a  record  obtained  when  a  rarefaction  moves  over  the 
foil  after  the  initial  shock  front  passage .  The  grid  lines  are  displaced 
back  towards  their  original  position,  along  the  line  CD,  indicating  a  de¬ 
crease  in  particle  velocity.  Reduction  of  such  foil  data  yields  the  relief 
wave  velocity  and  the  change  in  the  particle  velocity. 

Consider,  for  example,  the  measurement  of  the  relief  wave  veloc¬ 
ity  in  the  driver  using  this  technique.  Let  the  shock  in  the  driver  be 
produced  with  a  flying  plate  of  the  same  material.  The  shock  propagates 
into  the  driver  and  then  into  the  liquid  cell  where  foil  motion  is  observed. 
The  backward-facing  shock  produced  at  impact  is  reflected  at  the  rear 
of  the  flyer  plate  and  becomes  a  forward-facing  rarefaction  fan.  Figure 
12  is  a  time-distance  diagram  illustrating  the  wave  trajectories.  The 
thickness  of  the  cell  bottom  or  driver,  the  foil  angle,  and  depth  are 
chosen  so  that  the  head  of  the  rarefaction  wave  will  overtake  the  forward¬ 
facing  shock  front  when  it  is  still  interacting  with  the  foil.  Measuring 
the  depth  at  which  the  relief  wave  overtakes  the  shock  locates  the  point 
P  in  Fig.  12.  The  foil  analysis  gives  the  particle  velocity,  sound  speed, 
and  shock  velocity  in  the  fluid  so  that  segments  OP,  OS,  and  SP  can  be 
drawn.  The  points  A  and  B  can  be  located  with  a  knowledge  of  the  shock 
velocity  in  the  driver.  The  minimum  relief  wave  velocity  is  now  just 
the  slope  of  a  straight  line  segment  from  B  to  O.  The  absolute  magni¬ 
tude  can  be  computed  only  if  the  sound  velocity  in  the  driver  is  known  in 
the  triangular  region  COS.  By  using  thicker  ceil  bottoms,  the  error  in 
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FIG.  12  TIME-DISTANCE  DIAGRAM  ILLUSTRATING  THE  WAVE  TRAJECTORIES 
IN  SHOT  11,762 

the  initial  elastic  sound  speed  due  to  the  uncertainty  in  the  sound  veloci 
ity  along  segment  CO  can  be  reduced.  However,  since  the  elastic  relief 
wave  must  be  observed  in  the  fluid  cell  before  overtaking  the  shock  front, 
a  maximum  driver  thickness  exists. 

Several  problems  have  been  encountered  in  the  application  of 
the  immersed  foil  gage.  First,  the  use  of  an  evacuated  region  between 
the  flyer  and  driver  results  in  bowing  of  the  driver  somewhat,  depend¬ 
ing  on  its  construction.  The  result  is  that,  upon  flying  plate  impact,  a 
nonplanar  shock  is  produced  in  the  driver.  Since  the  analysis  is  one¬ 
dimensional,  such  a  three-dimensional  shockwave  is  undesirable. 
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Another  problem  exists  when  studying  attenuation  shocks  with  a 
gage  made  of  a  material  of  lower  shock  impedance  than  the  specimen  to 
which  the  gage  is  attached.  The  interaction  of  the  rarefaction  fan  that 
that  is  reflected  back  into  the  specimen  at  the  specimen-liquid  interface 
(CS  in  Fig.  12)  and  the  forward  facing  rarefactions  in  the  attenuating 
input  wave  may  cause  the  specimen  to  spall.  Therefore,  the  useful 
observation  time  of  the  initial  wave  in  the  fluid  cell  is  limited  by  the 
arrival  of  the  spall  signal. 

b.  Results 

Water  was  used  in  the  immersed  foil  gage  experiments  per¬ 
formed  on  this  project  because  of  the  amount  of  equation  of  state 
data  available.  The  gages  were  mounted  as  shown  in  Fig.  10  for  the 
flying  plate  experiments.  The  cell  bottom  was  either  annealed  or 
hardened  2024  aluminum.  The  flying  plate  assemblies  were  similar  to 
the  ones  used  for  the  low  velocity  attenuation  experiments. 

Two  experiments  were  conducted  where  the  annealed  aluminum 
cell  bottoms  were  0.080-inch  thick  and  the  aluminum  flying  plates  were 
also  0.  080-inch  thick.  The  plates  reached  a  velocity  of  0.  133  cm/jlsec. 
Figure  11  is  an  example  of  the  photographic  record  obtained.  The  shock 
wave  that  intersected  the  immersed  foil  in  this  shot  was  flat-topped,  i.  e. , 
the  reflected  relief  had  not  yet  overtaken  the  shock  front.  Th*»  locus  of 
trace  end-points,  line  AB,  would  be  a  straight  line  if  the  shock  velocity 
were  constant  and  the  wave  were  plane.  However,  due  to  the  bowing 
of  the  cell  bottom  when  the  flying  plate  chamber  was  evacuated,  the 
center  of  the  cell  base  was  impacted  by  the  flyer  before  the  edges  and 
a  curved  shock  front  propagated  into  the  water.  The  curved  line  from 
A  to  B  is  the  result  of  this  nonplanar  shock. 

Notice  that  the  traces  along  the  line  CD  began  to  move  toward 
their  original  position.  This  trace  motion  indicates  the  propagation  of 
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a  rarefaction  wave  through  the  water.  It  relieved  the  pressure  rapidly 
at  first  and  then  more  slowly.  The  initial  decrease  in  fluid  particle 
velocity  corresponded  to  the  arrival  of  the  relief  wave  originating  when 
the  left  -facing  shock  in  the  flyer  reflected  from  the  rear  of  the  flying 
plate. 

It  was  hoped  that  the  trace  records  after  jump  wou.'l  yield  more 
information  on  the  shock  profile  as  a  function  of  time.  However,  due  to 
the  uncertainties  caused  by  the  curved  shock  and  the  indefinite  motion 
of  the  traces,  the  additional  pulse  shape  information  has  not  yet  been 
deduced  from  the  data.  Calculations  based  on  this  experimental 
arrangement  indicate  that  the  cell  bottom  will  spall..  The  recording 
times  in  the  experiment  were  long  enough  to  observe  the  spall  signal 
but  it  has  not  been  identified. 

Shot  11,253  involved  the  use  of  a  water  gage  in  which  two  reflect¬ 
ing  Mylar  foils  were  mounted  one  above  the  other.  The  1/8-inch  flying 
plate  impacted  a  1  /8-inch  annealed  2024  aluminum  cell  bottom  at  0.  127 
cm/jisec.  The  additional  bottom  thickness  somewhat  reduced  the  curva¬ 
ture  of  the  shock  front.  The  record  also  indicates  the  motion  of  the  re¬ 
lief  waves  into  the  cell.  The  depth  at  which  the  relief  overtook  the 
shock  front  was  0.  313  inch  in  the  water. 

As  discussed  earlier,  a  measurement  of  the  relief  wave  velocity 
in  the  aluminum  cell  bottom  can  be  made  more  precisely  if  the  cell 
bottom  is  thicker  but  not  so  thick  that  the  relief  wave  BC  overtakes  the 
shock  front  in  the  aluminum  rather  than  in  the  water.  Shot  11,762  was 
conducted  in  order  to  measure  the  initial  relief  wave  velocity  in  2024- 
T351  aluminum.  Figure  12  is  the  corresponding  time-distance  diagram. 
The  1/8-inch  flying  plate  acquired  a  velocity  0.  128  *  0.02  cm/fjsec 
and  impacted  the  0.  507-inch-thick  aluminum  cell  bottom.  Calculations 
based  on  the  experimental  results  indicate  a  minimum  initial  relief 
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wave  velocity  at  110  kbars  in  the  aluminum  of  0.77  cm/fisec.  The  actual 
velocity  must  be  greater  than  this  because  the  sound  speed  in  the  region 
COS,  which  is  already  somewhat  relieved,  is  less  than  the  initial  re¬ 
lief  wave  velocity  in  the  fully  stressed  material. 

A  calculation  was  performed  using  the  constant  Poisson's  ratio 
elastoplastic  theory  to  determine  the  sound  speed  in  the  region  COS  of 
Fig.  12.  It  was  assumed  that  the  head  of  the  rplief  disturbance  moved 
through  a  uniform  region  along  OC  where  the  stress  was  35  kbars.  The 
resulting  sound  speed  was  0.  553  cm/fisec.  Employing  this  value  to  cal¬ 
culate  the  elastic  sound  speed  along  BC  gave  cej  =  0.81  cm/fisec,  which 
is  in  excellent  agreement  with  the  0.  80  cm/fisec  obtained  by  measuring 
the  free-surface  velocity  of  aluminum  samples  as  described  later. 


SECTION  IV 


EXPERIMENTAL  RESULTS 


1 .  Aluminum  (Low  Velocity  Flyer  Plates) 

Results  of  experiments  in  which  both  the  flyer  plates  and  the  tar¬ 
get  plates  were  as-received  2024-T351  aluminum  are  shown  in  Fig.  13, 
where  half  the  free-surface  velocity  is  shown  as  a  function  of  the  thick¬ 
ness  of  the  target.  The  data  are  plotted  in  this  way  to  facilitate  compari¬ 
son  with  calculated  results.  Target  thicknesses  are  given  in  multiples 


FIG.  13  PEAK  PARTICLE  VELOCITY  v».  SHOCK  TRAVEL  !«/*0)  FOR  2024-T351 
ALUMINUM  IMPACTED  AT  0.12S  e»/ps«c 


of  the  flyer  plate  thickness,  x  ,  which  for  the  experiments  reported 

0 

here  was  0.  125  inch.  Two  separate  experiments  gave  particle  veloc¬ 
ities  of  0. 0615  cm/usec  for  x/x  =  1,  i.e.,  the  free-surface  velocity 
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was  0.  123  cm/ji*ec  when  the  targets  were  0.  125-inch  thick.  This 
agrees  well  with  the  flyer  plate  velocity,  which  ranged  from  0. 124  to 
0.  126  cm/psec . 

Comparison  of  the  experimental  and  calculated  results  shows 
that  the  two  agree  fairly  well  as  to  the  point  at  which  attenuation  begins, 
approximately  5.  5  plate  thicknesses.  Closer  agreement  can  be  obtained 
by  using  a  smaller  increment  in  the  calculations.  For  the  results  pre¬ 
sented  here,  20  increments  were  used  for  each  plate  thickness  (0.  32  cm). 
The  curve  for  the  computed  results  shows  a  flat  portion  at  seven  to  nine 
plate  thicknesses.  The  experimental  results  show  no  tendency  for  the 
curve  to  become  flat.  This  continuous  decrease  of  the  particle  veloc¬ 
ity  in  the  experimental  results  forces  a  revision  uf  the  stress-strain 
relations  used  in  the  computational  model. 

Results  obtained  previous  to  this  project  for  2024- T4  annealed 
aluminum  are  given  in  Fig.  14  along  with  the  results  shown  in  Fig.  13. 
Results  of  Shot  10, 227  agree  well  with  the  new  data,  while  those  of 
Shot  10,  354  do  not.  It  is  probable  that  the  flying  plate  velocity  was 
somewhat  higher  in  Shot  10,  354  than  in  Shot  10,  227  and  the  plate  veloc¬ 
ity  experiments.  7 


However,  the  results  at  night  and  nine  plate  thicknesses  for 
Shot  10.  354  show  that  the  decay  rate  of  peak  pressure  with  distance  was 
almost  aero,  as  was  also  observed  in  earlier  experiments.9  Preference 
in  this  situation  is  given  to  the  new  data  which  give  nc  evidence  of  a 
plateau.  The  comparison  indicates  that  annealing  2024  aluminum  makes 
no  measurable  difference  in  the  target  thickness  at  which  attenuation 
commences.  The  difference  in  behavior  at  eight  and  nine  plate  thick¬ 
nesses  could  possibly  be  due  to  the  use  ot  annealed  aluminum  in  the 


FIG,  14  PEAK  PARTICLE  VELOCITY  SMOCK  TRAVEL  (*/«„)  FOR  2024-T35’ 

ALUMINUM  AND  2GM  ANNEALED  ALUMtbJW  MP  ACTED  AT  0,125  cm  «»»c 
($ho*»  10,227  <m4  10,354  win  t«pur*Ml  fw  #»  cnrwmbd 


earlier  study  o;  to  the  use  of  different  techniques  for  measuring  the 
free -surface  velocity. 

Results  of  attenuation  experiments  using  annealed  1060  aluminum 
are  given  in  Fig.  15  for  both  the  current  and  the  e*  shots.  If  the 
point  for  which  x/x^  *  1  and  the  particle  velocity  of  Q.  0685  is  ignored, 
the  average  particle  velocity  in  the  absence  of  attenuation  was  about 
0.0655  cm/psec  (using  only  the  ne**  dataj.  That  t#,  the  average  shun 
velocity  was  about  0.  127  c m /pse  The  velocity  of  the  1060  projectile 
plate  was  about  0.  126  e  0.  001  cm  :»sec  for  these  new  experiments 
Hence  the  higher  particle  velocity  tor  thin  1060  aluminum  targets  as 
compared  with  2024  targets  (Figs.  15  ar.d  U|  ir  at  leas!  partly  due  to 
the  fact  that  the  1060  flyer  plates  acquired  greater  velocity  than  did 
the  2024  flyer  plates.  Attenuation  sta-ts  at  about  six  plate  thicknesses 
in  1060  aluminum. 


r»&  15  PEA*.  PARTICLE  VELOCITY  v*.  SHOCK  TRAVEL  («/*„)  FOR  1060  ANNEALED 
ALiAUNUM  IMPACTED  AT  0.125  cm/ji**- 

The  results  of  the  new  study  of  1060  aluminum  shows  iomt  indi- 
cation  of  a  flat  region  at  seven  and  eight  plate  thicknesses  similar  to 
the  fiat  part  of  the  calculated  curve .  The  same  trend  is  evident  in  the 
results  of  the  earlier  Shot  10,  55 i,  which  had  *  higher  flying  plate 
velocity  than  Shot  10, 226. 

The  value  of  the  stress  relief  provided  by  the  elastic  wave, 
in  Section  II.  is  sot  obtainable  from  the  data  for  either  1060  or 
2024  aluminum  ic.  a  reliable  way.  As  shown  in  Section  71.  the  value  of 
is  necessary  for  a  determination  of  the  yield  strength  when  the 
elastopUstic  equation  of  state  .*  assumed.  Figure  14  gives  an  indication 
that  the  decrease  in  particle  velocity  may  be  about  0.  009  cm/psec,  cor¬ 
responding  to  a  drop  in  peak  shock  strength  of  25  kbars.  The  data  ler 


the  1060  aluminum  give  a  more  definite  indication  of  a  plateau  in  the 
plot  and  hence  a  more  reliable  value  ox  the  decrease  in  stress  (ae-CTf). 
The  value  appears  to  be  10  kbars. 

Table  III  gives  the  complete  results  for  all  attenuation  experi¬ 
ments  performed  during  the  current  program.  The  data  for  the  low 
velocity  plate  slap  experiments  on  aluminum  as  well  as  the  high  velocity 
experiments  with  Teflon  and  aluminum  are  presented. 

2 .  2024-T351  Aluminum  (High  Velocity  Flyer  Plates) 

Two  experiments  were  performed  with  2024-T351  aluminum  in 
the  as -received  condition  for  which  the  flyer  plate  velocity  was  about 
0.32  cm/|ise:.  This  velocity  corresponds  to  a  pressure  of  about  345 
kbars  in  aluminum.  Results  of  the  shots  are  given  in  Fig.  16  aiong 
with  the  calculated  curves  using  a  fluid-type  equation  of  state  and  the 
constant  Poisson's  ratio  elastoplastic  model.  The  average  shim  veloc  ¬ 
ity  in  Shot  11,  824  for  specimen  thicknesses  out  to  4.  1  plate  thicknesses 
is  0.  324  cm/^sec.  The  projectile  velocity  for  the  shot  was  0.  316  ± 
0.005  cm/(isec,  Hence  the  free- surface  velocity  appears  to  be  slightly 
greater  than  the  velocity  of  the  flyer  plate.  The  calculations  were 
done  using  a  velocity  of  0.33  cm/fisec.  Shot  11,861  shows  that  the 
particle  velocity  has  been  reduced  below  that  predicted  by  the  fluid 
model,  but  nwt  as  much  as  predicted  by  the  elastoplastic  model.  The 
attenuation  started  between  4.  1  and  5.2  plate  thicknesses,  implying  a 
sound  speed  of  about  0  93  cm/ysec  behind  the  shock  front.  The  data 
indicate  that  the  attenuation  of  stress  is  about  65  kbars,  and  there  is 
no  assurance  that  this  is  due  entirely  to  the  elastic  wave. 

3 .  Pressure  Dependence  of  Shear  Modulus  and  Yield  Stress  of 

2024-T351  Aluminum 

The  speed  of  the  head  of  the  elastic  relief  wave  car.  be  determined 
by  locating  the  point  at  which  it  overtakes  the  shock  front,  point  M  in 
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FLYING  PLATE  ATTENUATION  MEAS 
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PARTICLE  VELOCITY - cn/uwc 
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FIG.  16  PEAK  PARTICLE  VELOCITY  v*.  SHOCK  TRAVEL  (x/*0)  FOR  2024-T351 
ALUMINUM  IMPACTED  AT  0.33  cm/^«.c 


Fig.  2. 18  This  point  corresponds  to  the  target  thickness  at  which  the 
peak  particle  velocity  first  begins  to  drop,  as  shown  in  Figs,  13  through 
16.  There  are  other  methods  for  determining  the  relief  wave  velocity. 
Al'tshuler  et  al. 16  describe  a  "lateral- relaxation"  method  and  an 
"overtaking-relaxation"  method,  the  Utter  being  similar  to  the  tech¬ 
nique  described  above.  Still  another  is  based  on  the  analysis  of  the 
fluid  gage  records  (Section  111-3).  Both  the  overtake  and  fluid  gage 
methods  were  used  in  the  experiments  in  which  the  lower  velocity  plates 
impacted  2024  aluminum  targets.  Fluid  gages  were  not  used  with 
higher  velocity  plates  or  with  1060  aluminum  targets. 

From  the  measured  elastic  wave  speeds  and  associated  free- 
surface  velocities,  together  with  independent  Hugoniot  equation  of  state 
measurements,  values  of  the  elastic  shear  modulus,  G,  can  be  derived 
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as  shown  in  Section  II.  These  values  and  other  quantities  derived  from 
the  experiments  are  listed  in  Table  IV. 

The  eiastoplastic  model  as  currently  formulated  predicts  -*n  initial 
separation  of  the  elastic  and  plastic  relief  waves  with  a  relatively  con¬ 
stant  region  between  them.  This  separation  should  show  up  as  a  step  in 
the  decay  curve  as  indicated  by  the  theoretical  curves  of  Figs.  13 
through  16.  The  drop  in  particle  velocity  to  the  step  is  a  measure  of 
the  yield  strength  of  the  shocked  state.  The  data,  however  do  not 
clearly  show  such  a  step  except  possibly  for  1060  aluminum  (Fig.  15), 
and  one  shot  (No.  10,  354)  on  annealed  2024-T351  aluminum  (Fig.  14). 
Consequently,  values  for  the  yield  strength  [more  exactly  (<Te  -crf)] 
that  are  clearly  reliable  cannot  be  directly  obtained. 

It  should  be  noted,  however,  that  within  the  context  of  the  theoret¬ 
ical  model  assumed, all  parameters  except  the  yield  strength  are  deter¬ 
mined  directly  from  experiment  (including  G),  Consequently,  to  the 
extent  that  the  model  is  correct  an  effective  value  for  Y  can  be  deduced 
by  trial  and  error  by  requiring  agreement  between  the  calculated  and 
the  measured  decay  curves.  1 

Values  of  G  and  (Yg+Yj)  determined  from  the  experiments  are 
shown  as  functions  of  the  specific  volume  in  Fig.  17.  The  zero  stress 
value  of  G,  0.287  Mbar,  and  the  approximate  value  of  YQ ,  0.0025 
Mbar,  are  also  shown  in  *he  figure.  The  possible  errors  in  the  data  at 
the  high  stress  point  are  so  large  that  no  conclusive  inferences  can  be 
drawn  about  functional  relationships.  However,  several  relationships 
can  be  assumed  that  are  consistent  with  the  data;  these  permit  calcu¬ 
lations  of  shock  decay  to  be  made  that  can  be  usefully  compared  with 
experimental  attenuation  data. 

The  functions  chosen  for  an  initial  trial  comparison  of  theory  and 
experiment  are  the  curves  labeled  I  in  Fig.  17.  The  shear  modulus 
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Table  IV 

VALUES  OF  EXPERIMENTAL  PARAMETERS 
FOR  2024-T351  ALUMINUM 


Flyer 

Velocity 

Plate 

(cm/jj,sec) 

Parameters 

0.125 

0.33 

Peak  stress  (Mbars) 

Sound  speed,  c(cm/psec)* 

Sound  speed,  c(cm/psec)t 

0.110 

0.80  ±  0.02 
0.81 

0.345 

0.93  ±  0.05 

G  (Mbars) 

0.54  ±  0.07 

0.59  ±  0.25 

K  (Mbars) 

-  cjf  (Mbars) 

Ye  +  Yf  (Mbars) 

Coordinate  of  point  M,  Fig.  2 
Flyer  plate  thicknesses  (cm) 

1.27 

0.025 

0.013  ±  0.008 
5.5 

0.32(nominal) 

2.28 

0.065 

0.025  ±  0.0C8 
4.5 

0 . 32(nominal) 

*  Aluminum  free-surface  velocity  vs  depth  measurements, 
t  Immersed  foil  water  gage  measurement. 


(2 


(21) 


The  hydro stat  is  then  given  by 

P  =  0.764^  +  1.37m*  +  1.  103  m3 

The  curves  labeled  II  and  III  are  modifications  of  the  above  equa¬ 
tions  chosen  in  an  attempt  to  improve  the  fit  over  that  given  by  curves  I. 
These  are  primarily  alternative  functions  for  the  yield  strength,  since 
it  is  subject  to  the  largest  experimental  uncertainty.  A  weak  coupling 
between  the  shear  modulus  and  the  yield  strength  exists  for  a  given 
hydrostat;  this  coupling  is  the  source  of  the  differences  in  curves  1,  II, 
and  III  of  Fig.  17a. 

Note  that  the  functions  for  G  exhibit  a  maximum  value  near 
M  =  0.  217  (V  =  0.  295  cc/g).  The  corresponding  pressure  is  about  250 
kbar.  It  would  be  of  considerable  interest  to  better  determine  experi¬ 
mentally  whether  the  shear  modulus  possesses  such  a  maximum  along 
the  Hugoniot  curve,  since  this  would  indicate  a  trend  toward  true  fluid 
behavior. 

4.  Flow  Calculations  for  2024-T351  Aluminum 

Results  of  flow  calculations  using  the  assumptions  mentioned 
above  are  given  in  Fig.  18  for  the  lower  velocity  impact  case,  and  in 
Fig.  19  for  the  higher  velocity  case.  The  a,  b,  and  c  portions  of  each 
figure  refer  to  the  fits  designated  as  I,  U,  and  III  of  the  preceding 
section. 

For  lower  impact  velocities,  fit  I  (Fig.  18a)  shows  reasonable 
agreement  but  exhibits  a  stepwise  decrease  in  particle  velocity  that  is 
not  evident  in  the  data.  At  higher  velocities  (Fig.  19a)  this  fit  com¬ 
pares  favorably  at  five  plate  thicknesses  but  falls  off  too  quickly 
thereafter. 

Fit  II  shows  less  of  a  step  in  the  decay  curve  for  lower  velocity 
impact,  but  falls  off  too  slowly  at  the  greater  target  thicknesses  (Fig. 
18b).  The  agreement  for  higher  impact  velocity  is  quite  good  (Fig.  19b). 
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The  effect  of  increasing  the  number  of  cells  in  the  calculations 
is  shown  by  the  curve  labeled  "40"  in  Fig.  18b.  In  this  calculation  the 
flyer  plate  was  aoned  to  contain  40  cells  rather  than  20  as  used  in  all 
the  other  calculations.  It  may  be  noted  that  the  elastic  sound  speed 
(0.8  cm/fisec,  Table  IV)  was  determined  by  using  the  experimental 
value  of  5.5  plate  thicknesses  for  the  depth  at  which  *he  rarefaction 
overtakes  the  shock  front.  Zoning  the  flyer  plate  with  only  20  cells 
places  the  apparent  point  of  overtaking  at  x/*o  =  4.  5;  with  40  cells  this 
point  moves  to  x/xq  =  5.0.  Presumably,  convergence  to  the  value  5.  5 
would  occur  with  increasingly  fine  zoning. 

The  fit  Snown  as  curve  III  in  Fig.  17  clearly  gives  the  least  satis¬ 
factory  fit  to  the  decay  curves,  as  shown  in  Figs.  18c  and  19c. 

The  results  of  these  attempts  to  fit  the  decay  curves  indicate  that 
the  elastoplastic  theory  as  formulated  is  probably  oversimplified.  No 
step  in  the  decay  curves  can  be  clearly  identified,  at  least  for  2024-T351 
aluminum.  This  implies  that  there  is  no  pronounced  separation  of  the 
elastic  and  plastic  rarefaction  waves. 

The  most  likely  explanation  for  this  difference  is  that  a  Bauschinger 
effect  tends  to  spread  the  elastic  rarefaction  so  that  it  merges  with  the 
following  plastic  wave.  Bauschinger  effects  have  been  observed  in  plane- 
shock  waves  at  much  lower  pressures.  3'  * 

5.  Teflon 

The  attenuation  experiments  on  Teflon*  were  conducted  using 
thin  (0.0008-inch)  alumimsed  Teflon  shims  to  measure  the  initial  free- 
surface  velocity.  The  samples  were  arranged  in  the  name  manner  as 

*  Teflon  manufactured  by  AVCO  Corporation,  Wilmington.  Mass. 

Mean  density.  2.  195  g/cc. 
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the  aluminum  samples  described  in  Section  III  and  were  impacted  by  a 
1  '3 -inch  aluminum  flying  plate.  The  photographic  record  for  Shot 
11,864  is  shown  in  Fig.  20. 

The  four  Teflon  samples  were  of  different  thicknesses,  corres¬ 
ponding  to  multiples  of  flyer  plate  thicknesses  of  2.1,  4.1,  6.1,  and 
8.2.  In  Fig.  20,  the  line  AA  represents  the  time  at  which  the  shim 
first  moved.  Line  BB  indicates  the  shim  arrival  at  the  cover  glass 
after  having  crossed  a  0.48-cm  gap.  The  third  event  a*.  CC  is  inter¬ 
preted  as  the  arrival  of  the  surface  of  the  Teflon  sample  as  was  done 
earlier  for  aluminum  (Fig.  9).  The  presence  of  the  second  arrival 
indicates  that  the  Teflon  also  had  a  tenaile  strength  after  being  shocked. 

The  results  of  the  free- surface  velocity  measurements  appear  in 
Table  III.  The  flying  plate  velocity  was  0.  315  a  0.  010  cm/fisec  .  A 
graphical  impedance  solution  indicated  that  the  peak  induced  particle 
velocity  in  the  Teflon  was  0.  198  a  0.005  cm/fisec.  The  corresponding 

IT 

pressure  u-su  220  kbara. 


Note  that  thv"  thinnest  sample  was  only  two  flyer  plate  thickness. 
The  free-iuriace  velocity  measured  at  th  t  depth  was  0.439  a  0.  004 
cm/n»mc.  This  may  or  may  not  represent  the  peak  stress  in  the  Teflon, 
for  there  is  no  assurance  that  the  ehock  front  had  not  begun  to  attenu¬ 
ate  by  2x^  .  however,  a  time -distance  diagram  shows  that  the  rare¬ 
faction  wave  from  the  rear  of  the  flyer  plate  would  have  had  to  have  an 
almost  infinite  velocity  in  order  to  Overtake  the  shock  front  by  1^-x  . 

A  more  reasonable  rarefaction  yrave  velocity  predicts  the  initial  attenu¬ 
ation  to  occur  somewhat  deeper  than  2x  . 

0 


Considering  free- surface  velocity  of  the  thinnest  sample  to  repre¬ 
sent  the  amplitude  of  the  unattenuated  incident  shock,  it  is  seen  at  once 
that  the  free-surface  approximation,  i.e.  ,  u^  *  i“f,.  »•  •  poor  one  in 
this  experiment.  One  half  of  is  0.22  cm/paec,  whereas  the  induced 
particle  velocity  was  0.  198  cm/M»er 
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A  computer  calculation  was  performed  for  a  1/8-inch  aluminum 
plate  impacting  Teflon  at  0.  32  cm/jisec.  The  results  of  this  calculation 
are  shown  in  Fig.  21.  The  constant  Poisson's  ratio  equation  of  state 
for  aluminum  was  used  with  M  =  0.  055  (see  Section  II),  The  Teflon 
Hugoniot  was  represented  by 

PH  =  0.  112  4+  0.491  M8  +  0-248  <i3  .  (22) 

<L  8 

the  data  being  obtained  from  Netherwood. 


0*  M»«-40 


FIG.  21  CALCULATED  AND  EXPERIMENTAL  RESULTS  FOR  0.31-cm  ALUMINUM  PLATE 
IMPACTING  TEFLON  (Exporimantal  values  ara  1/2  uft  at  maaturad  in  Shot  11,864) 
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The  sound  speed  in  the  Teflon  was  then  obtained  from 

c  =  \j  dPH/dp  (23) 

where  Pj_j  is  the  pressure  along  the  Hugoniot. 

The  experimental  points  plotted  in  Fig.  21  are  one -half  the 
measured  free-surface  velocities.  These  show  strikingly  good  agree¬ 
ment  with  the  theoretical  curve,  but  to  some  extent  the  agreement  is 
fortuitous  because  the  free-surface  approximation  is  inaccurate. 
Nevertheless,  no  pronounced  elastoplastic  effects  are  evident  and  it 
appears  that  the  fluid  model  may  be  adequate.  More  information  on  the 
complete  P-V-E  equation  of  state  for  Teflon  would  permit  a  more  criti¬ 
cal  evaluation  of  the  fluid  model  to  be  made. 
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SECTION  V 
CONCLUSIONS 


Shock  wave  attenuation  data  are  presented  for  aluminum  of  types 
2024-T351  and  1060,  and  for  Teflon.  The  shock  waves  were  produced 
by  impact  of  aluminum  flyer  plates  accelerated  explosively  to  velocities 
of  0.  126  cm/jjsec  and  0.33  cm/jfsec.  Corresponding  peak  pressures 
in  aluminum  were  110  and  345  kbar,  respectively.  Free-surface 
velocities  were  measured  for  varying  target  thickness  by  observing 
with  a  streak  camera  the  time  of  flight  of  a  thin  shim  across  a  known 
gap.  This  technique  gave  higher  precision  than  techniques  employed 
previously. 6  » •> 7 

Reasonably  accurate  measurements  of  the  elastic  sound  speeds 
of  the  shocked  states  were  obtained  for  2024-T351  aluminum;  these  per¬ 
mit  values  of  the  elastic  shear  modulus  to  be  obtained. 

Attempts  to  measure  a  complete  particle-velocity  profile  by 
means  of  a  fluid  gage  were  not  successful.  Serious  complications  in 
the  analysis  of  the  experiments  were  introduced  by  changes  of  the  re¬ 
fractive  index  with  pressure  and  by  possible  spallation  near  the  speci¬ 
men  surface.  However,  an  independent  measurement  of  the  elastic 
sound  speed  was  obtained  at  a  pressure  oi  110  kbar  in  2024-T351 
aluminum  that  agrees  very  well  with  that  deduced  from  free-surface 
measurements. 

The  data  generated  on  2024-T351  aluminum  during  the  report 
period  do  not  show  a  stepwise  decrease  in  the  free-surface  velocity  as 
predicted  by  the  elastoplastic  model  with  a  von  Mines  (or  Coulomb) 
yield  criterion.  Instead,  the  decay  of  peak  velocity  is  continuous,  im¬ 
plying  that  there  is  no  distinct  separation  of  elastic  and  plastic  release 
waves.  Earlier  experiments  with  Manganin  wire  gages  and  optical 
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wedges  showed  similar  behavior  ’  '  but  the  evidence  was  not  regarded 
as  conclusive  at  that  time. 

The  lack  of  separation  of  elastic  and  plastic  relief  waves  is  un¬ 
fortunate  since  the  stress  relief  due  to  the  elastic  wave  is  a  measure  of 
the  yield  strength  under  shock  conditions,  but  it  cannot  be  clearly 
identified  in  the  data.  The  most  probable  explanation  is  that  a  Bausehin- 
ger  effect  tends  to  spread  the  elastic  relief  wave.  This  effect  has  been 
observed  at  lower  pressures.3 

The  case  for  annealed  1U60  aluminum  is  less  clear.  The  data 
appear  to  show  a  step  region,  but  more  data,  from  thicker  targets,  are 
needed  to  establish  its  behavior. 

An  important  result  is  that  a  significant  difference  is  observed  in 
the  behavior  of  2024-T351  and  1060  aluminum.  The  1060  is  clearly 
more  fluid-like  at  high  pressure,  Thus  a  correlation  is  observed  be¬ 
tween  the  zero  pressure  and  high  pressure  yield  strengths. 

The  experimental  results  for  Teflon  agree  reasonably  well  with 
predictions  assuming  fluid  behavior;  no  pronounced  elastoplastic  effects 
are  observed.  The  agreement  shown  is  evidently  partially  fortuitous, 
however,  since  an  assumption  Invoked  in  the  data  reduction  is  not  com¬ 
pletely  accurate.  *  More  information  on  the  complete  P-V-E  equation 
of  state  for  Teflon  would  permit  a  more  definitive  appraisal  of  the 
accuracy  of  the  fluid  model. 


*  The  computed  shock  particle  velocities  are  compared  with  one-half 
the  measured  freersurface  velocities.  This  procedure  is  valid 
only  where  entropy  can  be  neglected  in  the  equation  of  state  and 
material  rigidity  can  be  neglected.  For  Teflon  the  results  indicate 
that  this  "doubling  approximation"  is  not  accurate. 


56 


DEFINITION  OF  SYMBOLS 
IN  APPENDIX  I 


nj  Index  of  refraction  of  unshocked  liquid 
Index  of  refraction  of  shocked  liquid 
X  Displacement  of  grid  image  in  object  space 

g  Angle  through  which  the  foil  is  turned  by  shock 

a  Angle  between  the  shock  front  and  the  unturned  foil 

d  Distance  of  the  grid  from  the  glass  cover  on  the  fluid  cell 

s  Distance  from  the  glass  cover  to  the  Mylar  foil  in  the  fluid  cell 

Up  Particle  velocity  in  the  shocked  fluid 

qg  "Apparent  velocity" --the  velocity  of  the  point  of  intersection 
of  the  shock  and  foil 

Shock  velocity  in  the  fluid 

Pj  Initial  density  of  the  fluid 

p^  Shocked  density  of  the  fluid 

J  Displacement  of  grid  line  as  measured  on  the  film  reader 

(jump  in  image  space) 

R  Shot  reduction  factor  *  ratio  of  distance  on  the  film  reader  to 
9  distance  at  the  shot 

D  Effective  standoff  distance,  n^d  +  s 
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APPENDIX  I 

IMMERSED  FOIL  ANALYSIS 

The  cell  containing  the  inclined  aluminized  Mylar  foil  is 
mounted  directly  on  the  driving  material  in  which  a  shock  is  induced. 

The  fluid  is  added  to  immerse  the  foil  completely.  Small  holes  in  the 
foil  insure  that  the  hydrostatic  pressure  is  the  same  on  both  sides  of 
the  foil;  this  is  necessary  because  any  bowing  of  the  foil  changes  the 
reflecting  surface  from  plane  to  spherical,  fhvs  invalidating  the  follow  * 
ing  analysis.  The  light  source  and  grid  are  arranged  so  that  the  streak 
camera  views  the  images  of  the  grid  sources  in  the  immersed  foil 
(Fig.  22).  The  angle  between  the  normal  to  the  foil  and  the  optic  plane 
is  assumed  in  all  that  follows  to  be  so  near  zero  that  it  can  be  neglected. 


FIG.  27  SCHEMATIC  VIEW  OF  IMMERSED  FOIL  GAGE  (The  angle  b«t»w 
the  (heck  front  and  the  fail  i(  3,  end  the  fail  is  turned  through  on  angle  fi. 
nt  *  the  initial  index  ef  refraction  end  is  the  shacked  refraction  index.) 
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As  the  shock  propagates  into  the  liquid,  it  passes  over  the 
foil  and  imparts  to  it  a  velocity  in  the  direction  of  the  shock  velocity. 
Because  of  the  small  thickness  of  the  foil,  the  velocity  of  the  foil 
soon  becomes  equal  to  the  particle  velocity  of  the  shocked  liquid. 

The  effect  of  the  moving,  bent  foil  is  to  displace  toe  images  of  the 
grid  sources.  As  shown  in  Fig.  22,  the  new  images  are  displaced 
toward  the  portion  of  the  foil  previously  shocked.  The  displacement 
of  the  grid  image  along  the  direction  parallel  to  the  slit  is  called  the 
jump,  and  is  denoted  as  X  . 

Figure  1 1  is  a  streak  camera  record  obtained  with  the  use  of 
an  immersed  foil  cell.  The  shock  entered  the  gage  and  was  first 
detected  at  A  .  As  the  shock  moved  through  the  liquid,  its  intersection 
with  the  foil  was  indicated  by  the  successive  terminations  of  the  origi¬ 
nal  grid  traces  and  the  appearance  of  the  displaced  traces.  This  locus 
of  end  points  is  along  the  line  AB.  The  original  and  displaced  positions 
of  Line  11  are  denoted  in  the  figure.  This  line  displacement  as  recorded 
by  the  camera  is  proportional  to  the  jump,  X,  Also  note  that  the  dis¬ 
placed  grid  sources  are  nearer  to  the  driver  specimenthan  the  originals. 
There  is  also  an  interval  of  time  when  both  the  original  and  displaced 
images  are  photographed  by  the  camera.  This  occurs  because  the  dis¬ 
placed  image  can  be  observed  in  the  turned  foil  prior  to  the  arrival  of 
the  shock  at  the  original  position. 

Figure  23  is  a  schematic  of  the  shock  front,  turned 
Mylar  .  and  the  fh.'d  cell  cover.  Ti*c  ray  from  the  grid  source  is 
refracted  at  the  top  of  the  cell,  again  at  the  sh&ck  front,  and  is  then 
reflected  off  the  Mylar  foil  back  through  the  shock  to  emerge  parallel 
to  the  optic  plane.  The  distance  the  emerging  ray  is  displaced  from 
the  original  grid  position  is  the  jump  X  (the  camera  is  considered  to 
be  effectively  at  infinity).  Referring  to  the  figure,  it  is  evident  that 
the  jump  can  be  written  &s  / 

/ 

X  *  dtanf  +  Sjtan6+  s^l  tan  (a  -  c)  -  tan  (a  -  a)/  f24) 
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SHOCK  MONT-rOS! 
IWT£#SCCTtO*l 


FK5.  23  SCHEMATIC  Of  RELATION  BETWEEN  SHOCK  FRONT 
AM)  FLU©  CELL  COMPONENTS 
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Write  Snell  *  Law  or  each  refraction 

•  inf  =  n j  sin  6  , 

njSine  =  n^sinc  ,  (<>%) 

nj  «in a  =  aina 

and,  by  consideration  of  the  geometry,  we  can  write: 

6  =  a  -  e 

b  *  0  -  a  +  »  (26) 

c  =  a  -  0  -  b 


We  wiah  to  obtain  a  relationship  bet  een  the  turning  angle  ft  , 
the  shock  angle  a  ,  and  the  jump  X  .  Note  that  when  a  displaced  grid 
image  is  first  seen  in  the  turned  foil,  the  depth  of  the  shocked  fluid 
between  the  foil  and  the  shock  is  much  less  than  the  depth  of  unshocked 
flu  vj  between  the  shock  and  the  cell  cover  glass,  so  that  the  error 
introduced  by  replacing  s]  by  s  is  small.  It  is  also  assumed  that 
*.he  shocked  index  is  a  constant  in  the  region  between  the  shock  front 
and  the  turned  foil  and  that  can  hs  neglected.  Then  £q.  (24)  becomes 

X  ~  dtacf  +  s  tan  $  (<:7) 


Ueicg  Eqs,  (25)  and  <26) 


tVh'?:.  6  is  substituted  in  Eq.  128)  there  result*  the  desired  expression 
relating  X  ,  t»j  .  t>2  ,  4  .  and  6  . 
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The  shock  moving  through  the  fluid  imparts  a  particle  velocity, 

u  ,  to  the  fluid  which  is  parallel  to  the  shock  velocity,  U  .  Since 
P  s 

the  foil  moves  with  the  fluid,  Fig.  24  can  be  constructed  where  is 

the  apparent  velocity,  a  is  the  angle  between  the  shock  and  original 
foil,  and  &  is  the  angle  through  which  the  foil  turns.  From  the  geom¬ 
etry  of  the  figure  we  write 


FIG.  24  RELATIONSHIP  BETWEEN  SHOCK  VELOCITY  U,, 

FLUID  PARTICLE  VELOCITY  u  ,  APPARENT 
VELOCITY  q0,  SHOCK  ANGLE  a,  AND  TURNING  ANGLE  fi 


The  set  of  equations{28-3  1 )  gives  the  relationships  among  the 
parameters  that  describe  the  interaction  of  the  shock  with  the  liquid 
and  the  foil.  The  apparent  velocity  q^  and  the  jump  X  are  both 
obtained  from  an  analysis  of  the  experimental  film  record,  while  the 
grid  positions  and  cell  dimensions  are  measured  during  assembly.  'It 
is  then  Possible  to  determine  either  the  particle  velocity  or  the  re¬ 
fractive  index.  !f  a  knowledge  of  the  refractive  index  is  desired,  then 
the  fluid  particle  velocity  must  be  inferred  by  an  independent  measure¬ 
ment;  likewise,  when  using  the  technique  to  determine  particle  velocity, 
the  refractive  index  must  be  known  as  a  function  of  the  shock  state 


parameters. 


>mx:  ation 


When  an  enperinien*  is  performer,  in  which  the  angles  of  interest 
in  the  foil  analvsia  are  small,  i.  e.,  sin  6  =  tangent  0=0,  the  follow¬ 
ing  equations  are  applicable. 


X  ~  In^d  +  s)  5 


5  «  S7 

6  ‘  a  [‘ '  4] +  4  (trZrp) 


up =  V 


Combining  t2qs.  (32),  (33),  and  (34) 


/  nA  ni  qQX 
up  "  Us  +  n^  2  (s  4-njd) 
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It  has  been  found,8’  18  that  the  dependence  oi  refractive  index 
upon  density  is  very  nearly  linear  for  gome  liquids. 


n2  =  nl  +  k  *p2  "  °I* 


(36) 


From  one -dimensional  shock  theory 


«  pi 

t  *  '’~z 

Using  Eqs.  (36)  and  (37),  Eq.  (35)  becomes 


(37) 


u 


P 


=  QqX  /  <  2(n^d  f  s)  f  1 


VJ 


(38) 


Therefore  in  an  experiment  either  u  or  k  =  dn/dp  m a/  be  deter- 

P 

mined.  When  the  dependence  of  n^,  on  the  shock  parameters  is  not 
known,  then  the  Hugoniot  must  be  known  to  determine  x.^  .  Along  the 
Hugoniot,  u  =  u  (TJ  ),  and  Eq.  (35)  may  be  solved  for  n7  =  n,  (u  , 

p  p  8  w  c  p 

q^  ,  X  ,  etc. ).  The  particle  velocity  must  then  be  determined  by  an 
independent  means.  In  the  shots  where  the  flying  plate  velocity  was 
measured,  the  particle  velocity  in  the  aluminum  driver  was  known  and 
the  particle  velocity  in  the  fluid  was  found  by  impedance  mismatching. 
The  refractive  index  could  then  be  calculated  ior  che  initial  unatten¬ 
uated  shock.  The  dependence  of  the  refractive  index  on  density  was 
assumed  to  be  linear,  so  that  the  particle  velocity  could  then  be  inferred 
after  the  shock  had  begun  to  attenuate.  In  experiments  where  the  shocks 
were  attenuating  over  the  whole  region  of  measurement,  the  previously 
determined  refractive  index  function  was  employed. 

2.  Large  Angle  Analysis 

When  the  angles  associated  with  the  foil  cell  are  large,  so  that 
the  approximations  of  the  previous  section  cannot  be  used,  the  working 
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equations  are 


(28) 

(29) 

(30) 

(31) 

(39) 


(40) 

With  a  knowledge  of  refractive  index  and  Hugoniot  equation  of  state 
these  equations  can  be  solved  to  yield  the  particle  velocity  in  terms 
of  the  apparent  velocity  and  the  grid  line  jump. 

The  experiments  which  employed  large  foil  angles  were  designed 
so  as  to  measure  particle  velocity  over  a  substantial  depth  range  in 
water.  The  calibration  of  refractive  index  wa'  done  with  smaller  angle 
foils  and  the  results  were  applied  to  the  Urge  angle  experiments.  The 
relation 

n2  =  nj  +  k  (p2  -  Pj)  (41) 

was  assumed  to  be  sufficient.  Using  Eq.  (37),  Eq.  (41)  can  be  written 

n2  =  n^  +  k  Pj  («p/U§  -  tip)  (42) 


Now  define 


z  =  —  (arcsin!  —  since  I  4-  arcsin 

2I  n  ; 


SoEq.  (30)  becomes 


u  =  q  gfrJtt.-JS} 

p  ^0  cos  x 


X  =  n  jd  sin  6  |l  -  n  j  sin^  6  j 

h  r 

6  =  a -arcsint  ~  sin  2a- 

r*  L 


4-  s  tan  6 


2a  -  28  -  arcsin 


(?“■)]} 


Up  =  qQ  sin 8/cos  (a  -0) 
a  *  arcsin  (U  /qp) 
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Water  was  the  only  fluid  used  and  the  equation  of  state  was  that  of  Rice 
and  Walsh  80 


Ug  *  1.  483  +  25.  306  log  [  1  +  up/5.  190]  (43) 

where  U  and  u  have  units  of  mm/psec. 
s  p 

In  order  to  obtain  a  solution  to  the  set  of  equations  above,  a 

computer  program  was  developed.  The  program  solved  Eq.  (28)  by 

an  iterative  process  in  order  to  find  6 .  Then  by  choosing  a  value  for 

u  it  calculated  in  turn  U  ,  n„,  a,  and  z  and  compared  the  result 
P  s  2 

of  Eq.  (40)  with  the  initially  chosen  particle  velocity.  This  process 
was  repeated  until  the  percent  difference  between  the  chosen  and  cal¬ 
culated  particle  velocity  was  less  than  0.  1  percent. 

3.  Film  Record  Reduction 

Two  quantities  of  interest  are  obtained  from  the  photographic 
record:  (1)  the  grid  jump  and  (2)  the  velocity  of  the  shock-foil  inter¬ 
section,  i.e.,  the  apparent  velocity  (see  Fig.  11). 

The  jump  is  determined  by  measuring  the  distance  on  the  film 
that  the  line  is  displaced  and  dividing  by  the  appropriate  reduction 
factor  for  the  particular  camera  and  film  reading  device  used.  The 
films  on  this  project  were  read  on  a  Telereadex  instrument.  The  out¬ 
put  of  this  machine  is  in  counts  proportional  to  distances  on  the  film. 

It  is  then  necessary  to  measure  a  known  distance  on  the  shot  in  order 

to  determine  the  shot  reduction  factor  R  . 

s 

R  -  Counts  on  the  Telereadex 
s  ~  Distance  in  mm  on  the  shot 

The  distance  to  the  camera  objective  is  greater  for  the  grid  images  than 
for  the  shot  itself;  consequently,  the  reduction  factor  for  the  grid 
images,  Rg ,  is  somewhat  smaller  than  R#.  If  L  is  the  distance  from 
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D  the  effective  standoff  distance  for 


L  +  D 
— >*— 

The  magnitude  of  the  jump  as  measured  on  the  Telereadex  is  J . 

Therefore  the  actual  jump  magnitude  is  X  =  J/R  . 

6 

The  apparent  velocity  is  inferred  by  a  measurement  of  the 
slope  of  the  line  determined  by  the  ends  of  the  undisplaced  traces 
(see  Fig.  11).  In  general,  it  depends  on  the  coordinate  Z,  or  the 
distance  from  the  apex.  The  angle  between  the  vertical  and  the  tangent 
to  the  cutoff  curve  is  Y  .  Then  the  apparent  velocity  is 

-  q„  ••  jrcotY 

where  W  is  the  streak  camera  writing  rate  in  counts/psec  as  measured 
on  the  Telereadex. 

When  a  film  is  analyzed,  the  positions  of  the  undisplaced  and  dis- 
placed  traces  are  recorded.  These  sets  of  data  are  then  checked  by 
first  and  second  differences  and  smoothed  by  a  computer  program.  The 
program  also  differentiates  the  data  using  three  points  and  five  points. 

It  prints  out  the  cot  y  so  that  the  apparent  velocity  can  be  calculated  and 
the  procedures  for  determining  u^  as  described  above  can  be  carried 
out. 


the  shot  to  the  objective  lens,  and 
the  grid,  then  the  ratio  becomes 
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Appendix  II 

INDEX  OF  REFRACTION  MEASUREMENTS 
AND  DISCUSSION  OF  RELEVANT  DIELECTRIC  THEORY 

The  principal  goal  of  a  theory  of  dielectrics  is  to  relate  macro- 

» 

scopic  quantities  such  as  electric  polarization,  electric  field,  and 
electric  susceptibility  to  microscopic  concepts  and  quantities.  These 
relations  in  general  are  functions  of  temperature,  pressure,  density, 
and  the  frequency  of  the  applied  electric  field.  Much  work  has  been 
done  developing  both  formal  and  physical  theories  for  dielectric  media. 
The  purpose  of  this  appendix  is  to  review  the  basic  ideas  and  problems 
involved  in  order  to  appreciate  the  significance  of  new  data  that  have 
been  obtained  for  the  refractive  index  of  a  shocked  dielectric. 

1 .  General  Concepts 

A  quantitative  discussion  involves  the  so-called  electric  polarize- 

tion  P,  which  is  equal  to  the  electric  dipole  moment  per  unit  volume, 

and  its  relation  to  the  macroscopic  electric  field  E.  Provided  that  the 

medium  has  no  permanent  polarization,  that  is  that  in  the  absence  of  an 

external  electric  field  the  average  polarization  is  zero,  then  the 

polarization  is  a  function  of  the  applied  field  and  can  be  expanded 

as  a  power  series  in  the  field.  Experimentally  It  is  found  that  the 

linear  term  in  the  expansion  la  adequate  to  describe  physical  systems 

wherein  the  electric  field  is  not  comparable  to  interatomic  fields, 

0 

which  are  on  the  order  of  10  volts/cm.  If  the  further  simplification 
is  made  that  the  medium  is  isotropic  so  that  the  polarization  la  not 
dependent  on  the  direction  of  the  electric  field  with  respect  to  the 
dielectric,  we  can  write 


P  =  XpE  (44) 

where  the  proportionality  constant  is  called  the  electric  suscepti¬ 
bility.  The  first  problem  of  the  physical  theory  is  to  relate  the 
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macroscopic  quantity  xe»  or  its  equivalent  the  dielectric  constant  c,  to 
molecular  quantities.  Al tentatively,  since  x  i**  defined  bv  the  above 
relation,  the  theory  must  relate  P  and  E  to  molecular  quantities. 

The  calculation  of  the  local  field  is  not  simple;  in  fact,  one  of 
the  central  problems  is  to  relate  the  local  field  E '  to  the  macroscopic 
field  E.  In  only  a  few  special  cases  can  a  simple  relation  between  these 
two  quantities  be  derived.  Approximations  to  the  local  field  will  be  dis 
cussed  later. 

In  the  case  of  P,  the  problem  ia  simple.  Consider  a  molecule  in 
a  nonpolar  dielectric,  i.e.,  a  dielectric  in  which  the  molecules  do  not 
have  a  permanent  dipole  moment.  An  electric  field  exists  in  the  vicinity 
of  tne  molecule  due  to  an  externally  applied  field.  This  local  field  is 
denoted  as  E  to  keep  it  distinct  from  the  external  macroscopic  field  E. 
The  local  field  tends  to  polarize  the  molecule  and  the  induced  dipole 
moment  is  p.  Let  <pr>  be  the  mean  value  of  the  moment  computed  by  the 
methods  of  statistical  mechanics  so  that  the  desired  polarization  per 
unit  volume  la  just 

P  *  If  <p>  (4S) 

where  N  is  the  number  of  molebules  per  unit  volume. 

When  we  relate  microscopic  quantities,  we  still  cling  to  a  linear 
relation  between  the  local  field  and  the  extent  to  which  the  molecule 
is  polarized.  The  proportionality  constant,  a,  is  called  the  polariz¬ 
ability  and  is  a  microscopic  quantity.  The  molecules  may  have  a  perma¬ 
nent  dipole  moment  due  to  the  molecular  charge  distribution,  so  that 
the  total  moment  la  the  sum  of  the  permanent  moment,  JI,  and  the  induced 
moment,  a2'. 

p  *  JJ  ♦  ai f  ( 43 ) 

It  is  convenient  to  consider  a  aa  the  sum  resulting  from  the  contri¬ 
butions  of  several  sources  of  polarization  in  the  material.  The  main 
categories  are  labeled  electronic  polarization  or0,  and  atomic  polarlza- 

tion  a  ,  so  that  a  «  cr  +  a  .  A  concise  outline  of  the  several  sources 
"  a  e  a 

of  polarization  ia  given  by  Bottcher?*1 
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"The  average  (dipole)  moment,  caused  by  the  uniform  external  field, 
is  the  result  of  four  different  effects: 

a.  Rotation  (orientation)  effects. 

al:  The  uniform  (electric)  field  tends  to  direct  the 
permanent  dipoles  (dipole  orientation). 

a2:  It  also  tends  to  direct  an  anisotropic  particle 
into  such  a  position  that  its  axis  of  highest 
polarizability  coincides  with  the  direction  of 
the  external  field.  This  effect  may  be  neglected 
when  the  field  is  moderate. 

b.  Translation  (deformation)  effects. 

bl:  The  electrons  are  shifted  relative  to  the  posi¬ 
tive  charges  (electronic  polarization). 

b2:  Atoms  or  atom  groups  are  displaced  relative  to 
each  other  (atomic  polarization). 

The  rotation  effects  are  counteracted  by  the  thermal  movement  of 
the  molecules.  Thus  they  are  strongly  dependent  on  the  temperature, 
whereas  the  translation  effects  are  only  slightly  dependent  on  the 
temperature,  since  they  are  intramolecular  phenomena." 

These  classifications  are  convenient  because  the  dielectric  prop¬ 
erties  of  a  medium  are  dependent  upon  the  frequency  of  the  applied 
electric  field.  For  static  and  low  frequencies,  all  of  these  effects 
contribute  to  the  polarization.  As  the  frequency  is  increased  to,  ssy, 

the  microwave  region  for  water,  the  dipole  orientation  term  becomes 

9  ,  v 

very  Important,  At  about  2  x  10  cycles/sec  (X  =•  15  cm)  a  resonance 
occurs  which  corresponds  to  the  permanent  dipoles  of  the  water  flipping 
back  and  forth  in  phase  with  the  applied  field.  The  dipole  contribu¬ 
tion  becomes  less  Important  as  the  frequency  is  increased  further.  The 
molecules  require  a  characteristic  time  to  reorient  themselves  along 
the  applied  field,  and  when  the  frequency  of  the  applied  field  becomes 
too  large  it  changes  considerably  before  the  permanent  dipoles  are 
allowed  to  reach  an  equilibrium  distribution.  Hence,  at  optical  fre¬ 
quencies  the  effect  of  the  permanent  dipole  moment  on  tbe  polarization 
is  negligible. 
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Resonances  also  occur  for  the  atomic  effect.  Depending  on  the 
material,  resonant  absorption  due  to  atomic  oscillations  occurs  at 
infrared  or  optical  frequencies.  Similarly,  for  the  electronic  com¬ 
ponent,  resonances  appear  in  the  visible  and  ultraviolet.  Therefore, 
the  division  of  the  polarizability,  o,  into  its  components  Is  usefu' 
in  understanding  the  basic  physical  pbenome”*  that  determine  the 
dielectric  properties  within  a  particular  frequency  region. 

2.  Classical  Microscopic  Theory 

A  classical  treatment  of  the  problem  can  be  given  that  exposes 
some  of  the  essential  features.  Consider  the  displaceable  electrical 
charges  of  the  molecule  (these  charges  may  be  the  electrons  or  the 
atom  groups)  to  be  subject  to  a  restoring  force  that  la  proportional 
to  their  displacement 


Fj  -  -k trt  (47) 

The  electrostatic  forces  within  the  molecule  provide  the  restoring  forces. 
Let  so  electric  field  set  on  the  charges.  This  field  is  external  to  the 
molecule  but  is  Identified  with  the  local  field  in  which  the  molecule 

-e  g  th 

finds  itself,  E  .  The  equation  of  motion  for  the  1  particle  of 
charge  is 

d2r 

*,  *  V,  •  •>*'  ,4,) 

where  e^K'  la  the  external  force  on  the  1th  particle.  The  natural  fre¬ 
quency  of  the  bound  charge  is  so  that  the  equation  of  motion 

la  rewritten 


,  *  *  „  .  e  1 

1  ^  111  1 
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For  an  impressed  field  of  frequency  au/2rr  the  solution  is 


r  -  1  ?' 

1  ~  ,2  2  >  “ 
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- *) 


and  the  Induced  dipole  moment  is 


Pi  =  eiri 


51  ) 


so  the  polarizability  is 


,2  27 

BjCtDj  -  U,  ) 
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Since  the  polarizabilities  are  additive,  the  total  molecular  polariz¬ 
ability  is 


cr 


£ 

i 


~T~ 
■iK  ~ 


2v 

) 


'  S3 ' 


•■hen  the  frequency  of  the  external  field  approaches  one  of  the  natural 
frequencies  of  the  ays tea,  a  ■  u^,  a  damping  tern  aust  be  Included  in 
the  equation  of  notion.  This  treatment  indicates  the  Importance  of  the 
several  »'a  in  the  vicinity  of  a  resonance , 

This  classical  result  can  also  be  obtained  with  the  added  condition 
of  thermal  agitation.**  However,  an  exact  treatment  must  consider  other 
phenomena.  The  following  discussion  Is  by  V.  T.  Brown,  Jr.** 

"The  theory  of  dielectrics  would  be  simple  if  the  molecules. . . had 
no  electrical  effect  on  each  other  end  if  they  were  in  static 
equilibria.  In  feet,  however,  the  molecules  do  not  have  any  of 
these  characteristics.  There  are  three  characteristics  of  actual 
molecular  systems  that  complicate  the  problem:  (1)  they  obey 
quantise  lews:  (2)  they  interact  electrostatically:  and  ?3l  they 
are  subject  to  thermal  agitation. 

If  any  single  one  of  these  factors  is  introduced,  the  theory  be¬ 
comes  much  more  complicated.  Introduction  of  quantum  laws 
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(1)  alone  lead*  to  the  quentint-mechanlcal  formula  for  polariz¬ 
ability...;  to  evaluate  the  matrix  elements, . .rigorously ,  it  would 
be  necessary  to  solve  the  wave  'quation  for  each  type  of  molecule 
under  consideration,  and  this  is  a  formidable  problem  even  for 
only  moderately  complex  molecules.  Introduction  of  electrostatic 
interactions  (2)  alone  leads  to  the  problem  of  evaluating  the 
local  field  intensity  E7  lor  a  fixed  configuration  of  molecules; 
this  problem  has  been  aolved  rigorously  only  fci  crystal  lattices 
of  certain  types,  and  except  in  the  simplest  case  fa  cubic  crystal) 
the  calculation  requires  rather  elaborate  mathematical  techniques. 
Introduction  of  thermal  agitation  (3)  alone  leads,  however,  to 
a  compnratively  simple  problem,  one  that  can  be  solved  by  standard 
techniques  of  statistical  mechanics:  for  polar  molecules,  the 
solution  of  this  problem  gives  ue  Debye's  theory. 

If  any  two  of  the  complication  factors  are  Introduced  simultaneously, 
the  resulting  problem  c  Mid  eaailv  fill  a  book."21 «24 >26 

Even  though  the  situation  may  seem  hopeless,  we  will  show  that  with 
some  simplifying  assumptions  solutions  can  be  obtained  that  predict 
fairly  well  the  observed  phenomena. 

3.  The  Local  field 

We  have  shown  that  the  macroscopic  polarization  P  is  simply  related 
to  the  microscopic  dipole  moment  p  [see  Bq.  (45)].  The  local  field,  F.* 
is  going  to  be  the  resultant  of  several  fields.  First,  the  externally 
applied  field,  E,  contributes  to  the  local  field.  It  also  polarizes 
the  medium  so  that  there  will  be  a  contribution  to  the  local  field  of 
the  reference  molecule  due  to  cbe  polarization  of  its  neighbors.  If 
the  molecules  have  in  addition  a  permanent  dipole  aoaient,  the  local 
field  is  affected  by  dipole  Interactions  with  the  applied  field  sod 
with  the  dipole  field  of  the  reference  molecule.  *e  therefore  write 
the  local  field  as  the  eta  of  these  various  contributions. 


v  •  i  *  x  r 

i  1 

The  problem  la  to  evaluate  the  summation  using  some  reasonable  assump¬ 
tions  and  models.  The  macroscopic  relation  betmnmn  ?  and  1  [Bq.  '44)* 
la  recant  us*i«  the  definition  of  the  dielectric  constant. 


7* 


f  4  1  ! 
'54) 

(55) 

Equation  (55)  is  valid  for  isotropic  media  with  no  permanent  average 
polarization  (units  are  cgs). 

4 .  Low  Density  Model 

The  simplest  model  for  calculating  the  local  field  is  to  consider 
it  equal  to  the  external  field,  i.e.,  the  presence  of  the  polarized 
neighbors  has  a  negligible  effect  on  the  applied  field.  This  si tun  lion 
may  be  realized  in  a  low  density  gas.  Consider  the  molecules  to  have 
no  permanent  dipole  moment.  The  induced  dipole  moment  is  therefore 


Therefore 


p  - 

5  -  (1  ♦  ■lux,) 


— >  C*  —  1  — * 

p  J  ® - i  E 

4n 


<  p>  ~  a  E '  ~  a  E  v  ) 

The  polarization  per  unit  volume  is 

?  =  N  <y  E  ( 57 ) 

Utilizing  Eq.  (55)  the  molecular  poj arizability  is 


a  = 


g  -  1 
4ttN 


The  use  of  the  relation 


N 


(58) 


(59) 


where  is  Avogadro's  number,  p  is  the  mass  density,  and  M  is  the  molec¬ 
ular  weight,  results  in  ^he  expression 


( 4ttN  )  cv  =  -  (c  -  1) 
v  o  p  ' 


(60) 


I 
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In  order  to  give  meaning  to  Eq.  (60),  the  experimental  values  of  (e  -  IVp 

t 

might  be  examined.  If,  for  example,  this  quantity  is  constant,  one  may 
.justly  infer  that  a  is  also  a  constant  provided  the  approximation 
E*  =  E  is  valid  for  the  densities  involved.  1 

j 

This  model  can  be  extended  to  polar  molecules  with  the  use  of  statistical 

i 

mechanics.  We  have  for  polar  molecules 

i 

p  =  U  +•  0&'  (46) 

! 

The  usual  treatment  is  the  following,  First  the  contribution  to  the 
polarization  is  computed  as  if  p,  =  °>  which  gives 

e  -  1  =  (4t*)»  (58) 

Second,  the  contribution  dae  to  the  permanent  dipole  moment  is  calculated 
assuming  no  induced  dipole  moment  and  the  result  is 

t 

o  ’ 

e  -  1  =  (4nH)n  /3kT  I  (SI) 

Third,  these  two  contributions  are  added  to  give 


This  equation  is  called  the  Langevin-Debye  equation.  It  predicts  that 
for  polar  gases  the  quantity  (e  -  l)/p  should  be  linear  in  1/T.  This 
is  found  to  be  a  good  approximation  in  many  cases.  Notable  deviations 
are  observed  for  carbon  dioxide.28 

5.  High  Density  Model 

A  correction  to  the  previous  model  is  due  to  Lorentz,27  In  dense 
media,  the  local  field  is  no  longer  equal  to  the  applie^l  field  because 
the  polarization  of  neighboring  molecules  contributes  significantly  to 
the  local  field.  Consider  a  sphere  of  radius  "a"  about  the  molecule. 

The  distance  "a"  is  chosen  to  be  small  from  a  macroscopic  point  of  view, 
but  large  microscopically.  The  contribution  to  the  local  field  from 
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molecules  external  to  the  sphere,  B^,  can  be  calculated  by  using  macro¬ 
scopic  methods.  The  contribution  to  B'  from  molecules  within  the  sphere, 
Ug.  is  difficult  to  calculate,  for  there  we  are  forced  to  use  microscopic- 
methods. 

The  calculation  of  Bj  is  accomplished  by  considering  the  field  at  the 
center  of  a  hollow  sphere  of  radius  "a”  imbedded  in  a  dielectric  of 
uniform  polarization  P.  The  resulting  surface  charge  density  on  the 
sphere  surface  due  to  the  polarization  is  a  =  -P  •  n  where  n  is  i  normal 
to  the  surface.  The  contribution  to  the  local  field  at  the  sphere  center 
from  a  small  surface  charge  is 


d£ 


1 


(P  cos  8)  cos  0 

2 

a 


dS 


(63) 


where  the  direction  of  dE^  is  parallel  to  P.  Upon  integration  we  find 
that 


=  ~P  (64) 

The  total  local  field  becomes 

f'  =  B  +  |n  P  ♦  f2  (65) 

The  evaluation  of  B^  ^as  been  done  for  very  symmetric  situations  (Lorentz 
did  it  for  cubic  lattices)  and  found  to  be  zero.  When  the  molecules  are 
not  fixed  in  position,  Bg  is  not  generally  zero;  however,  for  this  model, 

the  approximation  is  that  B„  =  0.  Then  for  nonpolar  dielectrics,  the 

\  & 

local  field  becomes 


E'  =  B  +  P  (66) 

Using  Eq.  (55)  we  find  that  the  internal  field  is  related  to  the  external 
field  as 

V  («) 
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The  average  dipole  moment  is 


<p>  =  <*E '  =  c^-2— 5-^)  S  ( 68) 

The  polarization/unit  volume  becomes 

?  .  »<p>  .  1  (69) 

The  molecular  polarizability  is  obtained  by  combining  Eqs.  (68)  and  (69). 
The  result  is  - 

•  -  swyifH]  <-> 

where  N  is  defined  as  in  Eq.  (5S).  This  equation  is  known  as  the 
Clausius-Mosottl  equation. 

As  described  for  the  previous  model,  the  contribution  to  er  from  the 
permanent  dipole  moment  of  the  molecules  can  be  included.  The  result  is 

(*7«  5  •  ^  ♦  fe)  <"» 


6.  Optical  Frequency  Considerations 

In  order  to  relate  these  results  to  optical  frequencies,  the  Maxwell 
2 

relation  e  =  n  is  used.  This  equation  is  valid  so  long  as  the  magnetic 
permeability  of  the  dielectric  is  1.  This  can  be  seen  by  considering 
the  definition  of  refractive  index 

n  =  c/v  (72) 

where  v  is  the  phase  velocity  of  light  waves  with  frequency  f.  The 
speed  of  light  in  a  nonconducting  medium  is  c/^ii'c  where  is  the 
magnetic  permeability  of  the  medium.  Therefore 
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(73) 


Upon  substitution  of 


e  =  n  in  the  Clausius-Mosotti  equation  wc 


arrive  at  the  Lorentz-Lorenz  equation 


(74) 


Since  for  most  dielectrics  the  contribution  to  a  from  the  alignment  of 
permanent  dipoles  at  optical  frequencies  is  negligible,  we  are  justified 
to  drop  the  (j,  /3kT  term  in  Eq.  (71)  and  the  Lorentz-Lorenz  equation  can 
be  applied  equally  well  to  polar  molecules. 

A  quick  summary  of  the  validity  and  limitations  of  Eq.  (74)  is  given 
by  W.  F.  Brown.23 

"The  Lorentz  local-field  formula. . .has  been  derived  only  for 
a  lattice  of  dipoles  with  cubic  symmetry.  The  dipoles  must 
all  have  the  same  constant  vector  moment.  Actual  materials 
deviate  from  this  model  in  one  or  more  of  the  following 
respects:  the  molecules  are  not  located  at  points  of  a 

lattice;  the  fields  of  the  molecules  are  not  simple  dipole 
fields;  the  moments  are  not  all  the  same;  or  the  moments 
vary  (in  magnitude  or  direction  or  both)  with  time.  Con¬ 
sequently,  the  Lorentz  formula  and  the  dielectric-constant 
formula  that  follows  from  it... must  be  regarded  only  as 
first  approximations." 


Bottcher21  has  made  a  correction  to  the  Clausius-Mosotti  equation 
that  involves  the  "molecular  radius"  "a"  for  nonpolar  molecules.  His 
result  is 


9c 

2«  +  1)  - 


(75) 


Althoi«h  this  equation  predicts  the  behavior  of  some  gases  better  than 
the  Lorentz  equation,  some  calculations  by  Orttung2*  using  experimental 
data  indicated  that  there  was  not  a  clear  case  for  the  corrected  Bb'ttcher 
result  over  the  Lorentz  formulation  for  liquid  water. 


•• 
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The  molecular  refraction  la  defined  aa 


and  la  related  to  the  molecular  polarizability ,  a,  through  the  Lorentz- 
Lorenz  equation  (subject  of  course  to  the  limitations  of  the  Lorentz 
equation).  Since  we  are  dealing  with  optical  frequencies,  the  dominant 
component  in  the  polarizability  Is  the  electronic  component..  Examination 
of  R  as  a  function  of  density  for  a  particular  dielectric  may  Indicate 
the  extent  to  which  a  is  affected  by  changes  in  molecular  spacing.  To 
this  end,  Table  V  is  reprinted  from  the  Handbook  of  Physics.29  A  com¬ 
parison  of  R  is  made  for  large  changes  in  density  for  several  materials. 

Rg  and  R^  represent  the  molecular  refraction  for  the  gas  and  liquid 
phases,  respectively.  In  the  case  of  water  vapor,  R  =  3.74  cc/mole  and 

for  water  in  the  liquid  phase,  R  =  3.71  cc/mole.  This  corresponds  to 

X 

a  change  in  density  of  almost  100,000.  It  may  be  expected  that  an  addi¬ 
tional  increase  in  density  by  a  factor  of  less  than  two  in  a  shock  com¬ 
pression  will  not  affect  greatly  the  constancy  of  the  molecular  refraction. 

7.  Present  Experimental  Results 

Four  experiments  were  performed  (11,189;  11,190;  11,233;  and  11,762) 
in  which  an  aluminum  flying  plate  impacted  the  aluminum  cell  bottoms  of 
fluid  gages.  (See  Section  III  for  gage  description.)  The  flying  plate 
velocity  was  monitored  with  the  use  of  "pin"  switches.  On  Shot  11,762 

the  free  surface  velocity  of  the  aluminum  driver  was  also  measured  using 
the  shim  technique  described  in  Section  III.  The  peak  shock-induced 
particle  velocity  in  the  water  was  determined  graphically  with  the  use 
of  the  Rice  and  Valshae  Hugoniot  equation  of  state  for  the  water. 

Table  VI  lists  the  results  of  these  experiments.  The  index  of 
refraction  listed  and  the  shocked  state  densities  were  determined  from 
the  water  gage  records  with  the  use  of  the  computer  program  described  in 
Appendix  I.  The  dynamic  gage  measurements  were  all  completed  before  any 
rarefaction  waves  had  overtaken  the  shock  front.  The  records  themselves 
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Table  V 


MOLECULAR  REFRACTIVITIES  OF  VARIOUS  SUBSTANCES  CALCULATED 
FROM  GASEOUS"  AND  LIQUID-STATE  DATA* 


Substance 

B 

(ng  -1>  X  104 

D 

Pi 

Rg 

*1 

Hydrogen,  Hz 

2 

1.32 

1.10 

0.071 

1.98 

Water,  H20 

18 

2.49 

1 .334 

1 .000 

3.74 

Ammonia,  NH? 

17 

373 

1.325 

0.616 

5.60 

Nitrogen,  N2 

28 

2.96 

1 .205 

0.808 

4.34 

4.52 

Oxygen,  02 

32 

2.71 

1.221 

1.124 

4.06 

4.04 

Nitric  oxide,  NO 

30 

2.97 

1.330 

1.269 

4.45 

4.83 

Nitrous  oxide,  N.,0 

44 

5.16 

1.193 

0.870 

7.74 

6.24 

Chlorine,  Cl2 

71 

7.73 

1 .367 

1.33 

11 .58 

11 .92 

Hydrochloric,  HC1 

36.5 

4.47 

1.245 

0.95 

6.71 

5.95 

Bromine,  Bra 

160 

11.32 

1.659 

3.12 

17.00 

15.30 

Hydrobromic ,  HBr 
Sulfur: 

81 

5.73 

1.352 

1.630 

8.57 

10.72 

sa 

64 

11.11 

1.929 

2.04 

16.6 

14.9 

E,S 

34 

6.23 

1.384 

0.91 

9.35 

8.72 

so3 

64 

6.90 

1.410 

1.359 

10.32 

11.61 

csa 

76 

14.7 

1.628 

1.264 

22.0 

21 .4 

co. 

44 

4.49 

1.192 

0.796 

6.74 

6.78 

Methanol,  CHa0H 

32 

5.49 

1.331 

0.794 

8.24 

8.23 

Ethanol,  CaIfeOH 

46 

8.71 

1.3623 

0,800 

13.05 

12.72 

Acetaldehyde,  CR, CHO 

44 

8.11 

1.3316 

0.800 

12.16 

11.40 

Acetone,  CH,COCI^ 

58 

10.8 

1.3589 

0.791 

16.20 

16.05 

Phosphorus,  Pa 

62 

12.12 

2.144 

1.83 

18.20 

18.15 

*  The  data  refer  to  the  Ka  D  lines. 

Source:  Condon,  B.H.,  and  H.  Odlshav,  Handbook  of  Physics,  McGraw-Hill, 
Inc.,  New  York,  1958. 


State 

Refractive 

Index 

Molecular 

Refraction 

(cc/nole) 

ml 

Haul 

Density 

Cg/cc) 

Data  f roa  Present  Study 

Shot  nuaber 

11,762 

36 

175 

1.414 

1.493 

3.70 

11,253 

36 

175 

1.494 

3.88 

11,190 

39 

190 

1.437 

1.504 

3.71 

11,199 

38 

185 

1.433 

1.507 

3,76 

Atmospheric  Pressure  Data30 

Water 

Liquid 

1  X  10-* * * § 

20 

0.9962 

1.333 

3.712 

Ice 

1  X  io-* 

-3 

0.9164 

1 . 3090* 

3.777 

Vapor 

1  x  io-* 

0 

6  X  10-* 

1 .000249 

3.738 

Data  of  Ahrens  and  Rude  man® 

Shot  nuaber 

10,640 

.  «» 

240 

1.468 

1.474 

3.45 

10,363 

V  96 

395 

1.915 

1 .482 

3.39 

10,640$ 

49 

245 

1.4907 

1.5014 

3.4® 

Data  of  Zel'dovicb.  et  alJ 

19 

u»«e 

39 

190 

1.43 

1.47 

3.51 

110 

635 

1.67 

1.93 

3.33 

1 

i 

144 

660 

1.75 

1.56 

3.33 

•  H  *  16  ,01 5  g/aole 

t  TMpmtun  fraa  calculation  of  lie*  ud  Walsh,  Dmt .  90. 

t  Ordinary  my. 

§  tat*  reduction  by  praiut  author*. 


show  that  tha  initial  relief  waves  do  not  pass  over  the  Inserted  foil 
until  after  the  jumps  of  the  grid  traces  have  been  recorded. 

Table  VII  lists  the  molecular  refraction  of  water  as  calculated 
from  the  data  obtained  f row  the  experiments  mentioned  and  also  data 
obtained  from  two  other  sources.®*1*  The  average  value  of  R  for  our 
four  experiments  is 

R  3  3.71  ±  0.03  cc/mole 

ave 

at  a  density  of  1.43  ±  0.02  g/cc,  which  agrees  very  well  with  the  value 
of  R  computed  for  water  at  20*C  and  one  atmosphere. 

It  appears  that  the  electronic  component  of  the  aolecular  polariza¬ 
bility  for  water  is  little  affected  by  the  compacting  of  the  molecules 
froa  a  gaseous  state  to  the  densities  realised  at  38  kbars  in  the  shocked 
liquid. 

It  is  aeen  that  there  is  some  disagreement  between  our  data  and  that 
of  Ahrens  and  Rndenan*  and  Zel'dovlch  et  al.1*  The  disagreement  between 
our  data  and  that  of  Ahrens  nay  be  due  to  several  possible  sources: 
different  techniques  were  used  to  determine  the  shock- induced  particle 
velocity  la  the  water;  the  present  data  were  analyzed  using  the  Rice  and 
Walsh  equation  for  water  wfell*  the  earlier  data  ware  reduced  using  a 
Kugoalot  equation  derived  from  the  sane  experiment;  random  error. 

Zel'dovlch  et  el.  quote  a  ueaa  quadratic  error  le  V  of  ±0.01. 
However,  the  state  in  the  shocked  water  was  interred  by  a  asasuresent 
of  shock  velocity.  Small  errors  in  this  msssuremsat  might  account  for 
the  differences  between  their  data  sad  ours. 
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rapart  directly  bom  DOC  Otbor  fniifM  aaara 


(5)  “All  distribution  of  tbia  rapart  ia  contrail Qual¬ 
ified  DDC  aaara  ahali  matted  tbroonk 


If  tba  capart  baa  boon  (nalakaf  ic  tba  Office  af  Technical 
latitat.  Dapartmad  af  Cemmarce,  far  aala  ta  tba  pabltc.  lab 
cata  tbia  fart  ami  aatar  the  price.  if  larao 

It  SUPPLEMENT  ART  NOTES  Uaa  far  additional  sapiens. 


11  BOMMUSO  MOTT  ART  ACTTVITT:  Eotar  tba  anmr  af 
tba  fpdaaaHl  prajact  efftce  at  laboratory  ^caaonog  fyar 
b tf  far)  tba  taaaaacb  and  in  alapmaat  toclode  albaao 


IS.  ABB  TRACT  Eatar  an  abstract  giving  a  bwf  and  factual 
aaaamary  af  fa  dnenmarat  imdicativo  of  tba  rvpart.  avaa  though 
kt  aaap  alaa  app oar  alar  abara  ia  br  body  af  thr  irtkncrl 
pan.  R  afdtHraal  apaca  ka  itpliri.  a  cabunaattaa  abaat  aball 
ba  attatbad. 

B  ia  bMkty  daairabla  bn  tba  aba  tract  af  claaaifted  wprrta 
ba  traclaaaifkad.  Each  pataprapb  af  da  » bat  tact  aball  and  oritfe 
an  kndkraHrn  af  fa  military  a  act  tty  claaatfkcauao  af  tba  ta 
formats  an  is  tba  parogmgb.  tapraaanraf  aa  rr*».  f«i  rc>  »«  tVt 

Tbara  ia  aa  haitsuoa  co  tba  karmtb  af  tba  abstract  Ham- 
soar,  tba  usage  its  <  length  la  bam  1 30  to  333  a  ride 


co  tba  karmtb 
ka  bam  INtt 


14.  EET  WORD*  Eay  •orda  lit  tarbrri rally  nraatagfal  It  or 
or  abort  pbarrai  (bat  ckrracbtitr  a  (apart  aad  may  ba  aaad  t> 
indaa  aatttaa  far  crubpri  tba  (apart.  Eay  rratdr  mart  W 
aalarrad  aa  feat  a«  tacarity  cleaathcntlen  is  tapntad-  Idaati  - 
fkara.  aacb  aa  aaatpmaat  aaada!  dal ipnaWaa.  taada  aaara.  military 


prof coda  mama,  p  i  graphic  karat  tan.  may  ba  aaad  aa  bay 
•aria  bat  mill  ba  faUomnd  by  an  mdkrattan  af  (admiral  can 


toot,  tba  i 


af  tuba.  nbr.  aad  aadri  i»  optional 


aay  1km 


A rOC  nrr*  mat 


IscwRy  Cl—aiScnttaa 


